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ABSTRACT 

We present measurenients of the luminosity and mass functions of low-mass stars constructed from a 
catalog of matched Sloan Digital Sky Survey (SDSS) and 2 Micron All Sky Survey (2MASS) detections. 
This photometric catalog contains more than 25,000 matched SDSS and 2MASS point sources spanning 
~30 square degrees on the sky. We have obtained follow-up spectroscopy, complete to J=16, of more 
than 500 low mass dwarf candidates within a 1 square degree sub-sample, and thousands of additional 
dwarf candidates in the remaining 29 square degrees. This spectroscopic sample verifies that the 
photometric sample is complete, uncontaminated, and unbiased at the 99% level globally, and at the 
95% level in each color range. We use this sample to derive the luminosity and mass functions of low- 
mass stars over nearly a decade in mass (0.7 Mq > M* > 0.1 M©). The luminosity function of the 
Galactic disk is statistically consistent with that measured from volume complete samples in the solar 
neighborhood. We find that the logarithmically binned mass function is best fit with an Mc=0.29 log- 
normal distribution, with a 90% confidence interval of Mc=0.20-0.50. These 90% confidence intervals 
correspond to linearly binned mass functions peaking between 0.27 M© and 0.12 Mq, where the best 
fit MF turns over at 0.17 Mq. A power law fit to the entire mass range sampled here, however, returns 
a best fit of a=l.l (where the Salpeter slope is a = 2.35); a broken power law returns q;=2.04 at 
masses greater than log M — -0.5 (M=0.32 Mq), and a=0.2 at lower masses. These results agree well 
with most previous investigations, though differences in the analytic formalisms adopted to describe 
those mass functions, as well as the range over which the data are fit, can give the false impression 
of disagreement. Given the richness of modern-day astronomical datasets, we are entering the regime 
whereby stronger conclusions can be drawn by comparing the actual datapoints measured in different 
mass functions, rather than the results of analytic analyses that impose structure on the data a priori. 
Having validated this method to generate a low-mass luminosity function from matched SDSS/2MASS 
datasets, future studies will extend this technique to the entirety of the SDSS footprint. 
Subject headings: Surveys - stars: late-type - stars: low-mass - stars: luminosity function - stars: 
mass function - Galaxy: stellar content 



1. INTRODUCTION 

The mass function (MF) is a fundamental property of 
stellar systems, describing the number of stars as a func- 
tion of stellar mass. The MF is thus a statistical measure 
of the end result of the star formation process. Succinctly 
characterizing a stellar population, the MF also informs 
our understanding of the structure and dynamical evolu- 
tion of stellar clusters, the Milky Way and other galaxies. 

The fi r st me asurement of the MF was conducted by 
ISalpeten (|1955| ): for higher mass stars, this seminal re- 
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suit remains essentially unchanged to the present day. 
Salpeter found that the MF can be fit as a power law, 
formally expressed as: 



^(M) = ^ocM- 



stars pc 



Mr 



(1) 



where a = 2.35 is known as the 'Salpeter slope'. Power 
law MFs which find a larger fraction of low-mass stars 
have larger values of alpha and are expressed as being 
'steeper', while mass functions that find a lower fraction 
of low-mass stars have smaller values of a and are de- 
scribed as 'flatter'. 

Other functional forms have been suggested to charac- 
terize the MF. In particular, a num ber of investigators 
(jMiller fc Scalolll979l : IChabrieH[2005[) present log normal 
MFs, which can be expressed as 
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Whichever functional form they adopt, most studies 
agree that the stellar MF appears to reach a peak at a 
few tenths of a solar mass. Power law characterizations 
of the MF require a flattening from a Salpeter slope for 
M > 1 Mq to a ~0.5-1.7 below 1 Mq, while iChabneil 
()2005f ) finds a log-normal MF with a characteristic mass 
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of Mc=0.2 is required to reproduce the stellar MF for 
a volume complete sample in the solar neighbor hood. 
We re f er the reader to compreh ensive r eview s by IScalo 
(|T98 ^ . [Reid"fc HawlevI (2000). Krouoal (poM l. lChabrler 
(2003), and lCorbelli et a l. (2005) for a more detailed un- 
derstanding of the vast literature devoted to measuring 
the stellar MF. 

Several processes have been advanced to explain the 
shape of th e MF; examp l es inc lude gravitational frag- 
mentation ([Klessen et al] 11998) , competitive accretion 
(|Larsonlll992f ). truncation of m ass accret ion due to ra- 
diative or dynamical feedback (|Silk|[l995D . star-star in- 
teractions (jReipurth fc Clarkil2001[ ) and the primordial 
distribution of clump ma sses within molecular clouds 
(jPadoan fc Nordlundll200 2) . The efficiency of each mech- 
anism is influenced by physical variables, such as the 
metallicity and magnetic field strength of the parent 
molecular cloud, the local stellar density, and the inten- 
sity of the surrounding radiation field. These effects may 
ultimately result in observable MF variations as a func- 
tion of environment, but at present the observed varia- 
tions can be largely attributed to uncertainties arising 
from finite s ample sizes and systematic differences in 
methodology ()Chabrieill2005f) . though intriguing signs of 
MF variations with environment may be emerging (e.g., 
Taurus; Luhman 2004). 

Previous studies of the Galactic disk MF have been 
limited, however, to samples of a few thousand low-mass 
stars (with low-mass stars defined for the purposes of 
this paper as M* < 0.7 Mq), due to the inability to si- 
mult aneously obta i n a d eep and wide photometric sam- 
ple. iHawlev et al.l (|2002D demonstrated that Sloan Dig- 
ital Sky Survey (SDSS) and Two Micron All Sky Sur- 
vey (2MASS) photometry of low-mass stars shows mono- 
tonic behavior across a wide range of colors. The accu- 
rate, multi-color catalogs produced by SDSS and 2MASS 
thereby allow the identification and characterization of 
millions of low-mass stars in the local Galactic neigh- 
borhood, enlarging previous photometric samples of field 
low- mass dwarfs by several orders of magnitude. The sta- 
tistical power of such a sample makes an SDSS/2MASS 
catalog of low-mass dwarfs a promising avenue for im- 
proving our measurement of the low-mass stellar lumi- 
nosity and mass functions. 

To test and calibrate this technique before applying 
it to the entirety of the SDSS and 2MASS databases, 
we have performed a combined photometric and spec- 
troscopic study of nearly 30,000 stars detected in ~30 
square degrees of overlap between the SDSS and 2MASS 
footprints. We describe in ^ the construction of this 
sample by matching the SDSS, 2MASS, and Guide Star 
Catalog (GSC) photometric databases. In Sj3]we present 
the spectroscopic follow-up of more than a third of the 
sample. We verify the completeness limit of the photo- 
metric sample in §4, where we also analyze the spectro- 
scopic sample to understand the level of contamination 
and bias within the photometric catalog. We use this 
photometric catalog to measure the luminosity function 
of low-mass stars in the Galactic neighborhood in §5, 
derive a measurement of the low-mass stellar mass func- 
tion in §6, and summarize our conclusions in §7. An ap- 
pendix to this work presents native or transformed SDSS 
photometry for stars with measured trigonometric paral- 
laxes, useful for constraining the empirical SDSS/2MASS 



color-magnitude relation applied in §5. 

We note that a full and complete description of our 
analysis necessarily results in a rather lengthy paper; we 
therefore suggest that many readers may wish to skip 
directly to §5 fc 6 to concentrate upon the core scientific 
findings of this work, and then return to read §2, 3 & 4 if a 
deeper understanding of the assembly of the photometric 
sample and the effects of incompleteness, contamination 
and bias is desired. 

2. A UNIFIED CATALOG OF SURVEY PHOTOMETRY 

We analyzed SDSS, 2MASS, and GSC photometry for 
objects lying within a 'calibration region' defined as the 
area between right ascensions (a) of 3 and 15 degrees and 
declinations (S) of -1.25 to -f 1.25 degrees. An abnormally 
large number of SDSS spectra of point sources were ob- 
tained in this region, allowing a test of the efficiency and 
robustness of the photometric sample assembled here (see 
i )3.1[) . We excluded a small area [13.1 < a < 13.3, S < 
-1.05] surrounding HD5112, a bright (V = 4.77) MO gi- 
ant, which badly saturates SDSS photometry over a large 
area. The area of the calibration region thus subtends 
29.957 square degrees on the sky. 

2.1. Individual Survey Datasets 

2.1.1. Selecting 2MASS Stars 

The Two Micron All Sky Survey (2MASS; 
ISkrutskie et al.l l2006f ) provides a homogeneous cat- 
alog in three near infrared filters [JHKs) ranging from 
1 to 2.4 ^m. We queried the 2MASS All Sky Data 
Release (|Cutri et al.ll2003[ ) via the GATOR interface^ 
for all objects detected within the calibration region, 
retaining only objects with: 

• astrometric location within the boundaries of the 
calibration region; 

• high quality J band photometry (jphqual flag = 
'A' or jrdflag = '1' or '3'); 

• unique detections of catalog objects (use-src flag = 

'1'); 

• point source morphology (ext-key = 'null'); 

• no association with the predicted location of a 
known minor planet, comet, planet or planetary 
satellite (mp_flg = 0). 

These selection criteria generate a sample of 30,449 
stellar candidates with IR detections within the calibra- 
tion region. 

2.1.2. Selecting SDSS Stars 

The initial mission of the SDSS was to map a quarter 
of the night sky centered on the North Galactic cap, ac- 
quirin g accurate photometry of 100 million o bjects in 5 
filters (jFukugita et al.llT996l; iGunn et al.lll998l) and accu- 
mulating over 1 million spectra ( York et al.ll200Cl ). This 
mission is now complete, with over 9000 square degrees 
of photometry and over 1.25 million spectra released to 
the public in Data Release 6 (jAdelman-McCarthv et aLl 

^ http://irsa.ipac.caltcch.edu/applications/Gator/ 
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I2OO8D . SDSS observations were obtained with a dedi- 
cated 2 5 meter telescope at the Apache Point Obser- 
vatory (|Gunn et al.1 120061 ). Photometric data were ac- 
quired as the telescope's 3 degree field imaged in 5 fil- 
ters {u,g,r,i,z) nearly simultaneously by scanning a great 
circle across the sky at approximately the sidereal rate. 
Photometric data were reduce d by a custom SPSS data 
processing pipehne (PHOTO ; iLupton et al.ll200ll ) with 
calibrations obtained from observation s by a 20-inch pho- 
tometric telescope at the same site jHogg et al.l 120011 : 
iSmith et all l2002t iTucker et all 120061). The article de- 
scribing the Early Data Release ( jStoughton et al.ir2002f l 
provides information on the central wavelengths and 
widths of the SDSS filters, while papers discu ssing the 
"asinh" magnitude system (Lupton ct al." 1999) and the 
SDSS standard star system (S mith ct al. 200 2) provide 
further information on the calibration of the SDSS pho- 
tometric system. Astrometric pre cision and d ata qu ality 
assurance a r e also described by iPier et all (j2003[ ) and 
llvezic et al.] (|2004[ ). 

We assembled a catalog of stellar SDSS sources within 
the calibration region from imaging data firs t made pub- 
he as Data Release 2 (jAbazajian et al.l[200l ). SDSS stel- 
lar candidates satisfied the following criteria: 

• SDSS astrometric location within 0.005 degrees 
(15-20 ") of the calibration region (to allow matches 
for 2MASS objects whose proper motion may have 
carried them slightly out of the calibration region 
by the time of the SDSS imaging epoch) ; 

• observed during a run which passed basic data 
quality tests (GOOD flag set); 

• object represents the primary detection in the 
SDSS photometric database (PRIMARY flag set) 
to prevent multiple detections from overlap areas 
between adjacent SDSS runs; 

• morphological identification as a point source 
(TYPE = 6). 

• i and z magnitudes above survey completeness lim- 
its {i < 21.1, z < 20.3); 

• no photometric processing flaws to ensure ac- 
curate photometry (SATURATED flag not set 
in either r, i or z, PEAKCENTER flag not 
set, NOTCHECKED flag not set, DEBLEND- 
NOPEAK flag not set for an object with i psf er- 
ror > 0.2 magnitudes, PSF-FLUX-INTERP flag 
not set, BAD-COUNTS-ERROR flag not set, 
INTERP-CENTER and CR flags not both set); 

These selection criteria identified 76,966 high quality 
SDSS detections of stellar candidates within this region; 
the last set of quality cuts removed 1486 stellar candi- 
dates, or ~1.9% of the initial stellar catalog. 

We note a subtle source of uncertainty in the SDSS z 
band photometry which may be of more importance for 
this study than for most other uses of the SDSS database. 
Spectral synthesis of L and T dwarfs (Burgasser, private 
communication) indicates that slight differences in the 
spectral response of the CCD detectors used for the z 
band observations can introduce z magnitude variations 



up to 0.1 mags for early L through late T dwarfs. As we 
lack a robust characterization of this effect, and the vast 
majority of our sample is composed of K and M dwarfs 
for which this source of uncertainty is less important, we 
have not applied any z-band corrections to the survey 
photometry. 

2.1.3. Selecting GSC Stars 

The depth of the SDSS comes at the expense of a rel- 
atively faint saturation limit of isat ^ 14. This fixed ap- 
parent magnitude limit imposes a lower limit on the dis- 
tance to SDSS stars with good photometric detections, 
and this distance limit is itself a function of the abso- 
lute magnitude of the star. Typical MO dwarfs within 
100 parsecs are saturated in SDSS imaging, for example, 
while M8 dwarfs only saturate in i for distances less than 
10 parsecs. The 2MASS faint limit, however, restricted 
our sample to objects within 100 parsecs for objects of 
type MS and later. Constructing a luminosity function 
of stars detected in both SDSS and 2MASS then required 
that early and late M dwarfs be sampled from disjoint 
physical volumes. As the calibration region is located at 
a Galactic latitude of -62 degrees, this could potentially 
have introduced a bias in our analysis related to the scale 
height of the Galactic disk. 

To include stars with saturated SDSS imag- 
ing in our catalog, we made use of GSC 
(version 2.2) optical photometry. The GSC 
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(I Space Telescope Science Institute fc Osservatorio Astronomico di 1 
|2001| ) is based on scans of photographic plates from 
the Palomar and UK Schmidt telescopes and provides 
star/galaxy morphological classification from measure- 
ments of 500 million objects in photographic Bj (as 
defined bv lReid et al] 1 19911 A - 4500 A) and Rp (A ~ 
6500 A) to a completeness limit of Rp ~ 18.5. 

We assembled a catalog of GSC (version 2.2) objects 
brighter than a magnitude of 19.5 using the GSC data 
access page^'^^, which we searched for morphological point 
sources (GSC Classification code = 0). These criteria 
identified 33,612 stellar candidates with GSC detections 
within the calibration region. 

2.2. Creating a Matched Sample 

To generate a comprehensive optical/near-infrared cat- 
alog of calibration region stars, we have merged the 
2MASS, SDSS, and GSC steUar samples into a single 
matched catalog. As 2MASS photometry provides our 
only source of infrared observations, only objects with 
2MASS detections were retained. 2MASS stars were 
matched to SDSS and GSC counterparts using a maxi- 
mum matching radius of 5". We required unique matches 
- if an optical source was a potential counterpart for 
multiple NIR sources, only the closest association was 
preserved. 

Of the 30,499 2MASS sources in the calibration region, 
19,617 possessed both SDSS and GSC stellar counter- 
parts. An additional 5265 and 4608 2MASS detections 
were identified with only a GSC or SDSS stellar coun- 
terpart, respectively. The final 959 2MASS stellar can- 
didates (^ 3% of the input 2MASS sample) were not 
matched to counterparts in either optical catalog. Of 

^^^ http://www-gsss.stsci.edu/support/data_access.htm 
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these unmatched sources, 606 have 15.2 < J < 16.3, and 
753 have J < 16.3; given the J = 16.2 completeness hmit 
of our catalog (see M.l^ . this indicates that the bulk of 
these unmatched sources lie at the faint end of our sam- 
ple, but still are confident detections. Visual analysis 
of SDSS imaging of the 753 J < 16.3 sources identified 
73% (550/753) as galaxies: these objects were unresolved 
by 2MASS imaging, but resolved by the optical surveys, 
and were therefore not contained in the SDSS and GSC 
stellar catalogs. 

The remaining 203 unmatched 2MASS sources brighter 
than J = 16.3 were stellar sources, dominated by sources 
too faint for the GSC catalog and with sub-optimal SDSS 
photometry, and therefore excluded by the quality cuts 
described in ii2.1.2l As these unmatched stellar sources 
represent a possible source of incompleteness in our final 
sample, we explore their properties in more detail in ^A.l\ 

2.3. Recalibrating GSC Photometry 

The 19,617 sources in common between the SDSS and 
GSC catalogs allowed an assessment of zero point errors 
in the GSC data. We calculated synthetic GSC magni- 
tu des for SPSS sourc es using equations originally derived 
bv lSesar et al.l ()2006[ ). modified to account for local zero- 
point offsets: 

Bj,syn = .9 + 0.279(g - r) + 0.06 (3) 

RF,syn = r~ 0.209(5 - r) - 0.09 (4) 

The shape of the calibration region, 12 degrees wide in 
right ascension and only 2.5 degrees across in declination, 
implies that spatial gradients in imaging quality in right 
ascension are more important than gradients in declina- 
tion, an expectation borne out in our dataset. The top 
panel of Figure [T] shows residuals between observed and 
synthetic GSC magnitudes as a function of right ascen- 
sion. Assuming that the digital, drift-scanned SDSS data 
has a uniform zero point across the calibration region, 
these residuals reflect changes in photographic sensitiv- 
ity across POSS plates. This expectation is supported 
by sharp discontinuities in the residuals on 5 degree spa- 
tial scales, the size of the POSS plates from which GSC 
magnitudes are measured. 

To correct for spatial variations in the GSC zero point, 
we calculated the mean offset between synthetic and ob- 
served GSC magnitudes in 0.25 degree bins of right as- 
cension. Applying this offset to all objects in that spatial 
bin produces recalibrated Bj and Rp magnitudes, cor- 
rected for the presence of the known zero-point error. 
Though these systematic zeropoint offsets arc largely a 
function of right ascension, we have also repeated this 
process using 0.1 degrees bins in declination. The result 
of this recalibration is illustrated in the bottom panel of 
Figure [U which shows that large scale spatial gradients 
in the residuals have been removed. 

We adopt characteristic errors for the recalibrated GSC 
magnitudes based on the remaining residuals between 
the recalibrated and synthetic GSC magnitudes. Gaus- 
sian errors with a =0.13 accurately describe the recali- 
brated Bj residuals across the entirety of the calibration 
region. The recalibrated Rp residuals are well described 
by a Gaussian distribution with cr = 0.1 magnitudes for 
a > 8.2; at a < 8.2, a ^^0.15 is required to fit the bulk 
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Fig. 1. — Residuals of observed and synthetic GSC Bj mag- 
nitudes calculated from SDSS photometry as a function of a b 
efore (top) and after (bottom) recalibration with SDSS-based zero- 
points. The structure on 5 degree spatial scales in the top panel 
suggests that GSC magnitudes, measured fr om digitized POSS 
plates 5 degrees in size, contain systematic zero-point errors of the 
same order as the random photometric errors. 

of the sources, and a tail of sources with large errors 
is still present. We thus conservatively assign errors of 
0.2 magnitudes to recalibrated Bj magnitudes, and 0.15 
magnitudes to recalibrated Rp magnitudes. Given the 
increased uncertainty of the Rp magnitudes for a < 8.2 
degrees, we rely largely on Bj magnitudes in conducting 
our analysis. 

2.4. Synthetic SDSS/2MASS Colors for GSC Sources 

Measuring the luminosity and mass functions from 
photometric observations of the Galactic field requires 
the adoption of a color-magnitude relation to estimate 
the absolute magnitude, and ultimately mass, of each 
star. We used i — J for this purpose, as iHawlev et al.l 
(|2002f ) have demonstrated that this color is a monotonic 
indicator of spectral type, and therefore temperature and 
luminosity, for late type dwarfs. 

To create a sample of uniform i — J colors for the full 
photometric sample, we derived synthetic i~ J colors for 
objects with observed Bj — J colors. Figure [5] shows the 
Bj — J and i~J colors for sources detected in both optical 
surveys. Interpolating the Bj — J colors of our sample 
onto the median trend of i — J as a function of Bj — J, we 
assigned each GSC/2MASS object that lacked bona fide 
SDSS data a synthetic i — J color. We estimated errors 
by combining in quadrature photometric uncertainties in 
Bj — J and the dispersion of the Bj — J vs. i — J relation. 

Using the same methodology, we also derived synthetic 
i — z and r — i colors for all GSC sources with Bj — J 
colors. For GSC sources with only Rp magnitudes, we 
calculated the full suite of synthetic SDSS colors {i — J, 
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Fig. 2. — Bj — J and i — J colors for objects observed in both 
optical catalogs shown as black points and contours. The median 
i — J colors measured for 0.1 magnitude bins of Bj — J color are 
shown as triangles; the dashed line shows the function used to 
assign synthetic i — J colors as a function of Bj — J color to sources 
that lack SDSS data (e.g., saturated). 

r — i, i ^ z) using the same technique as appHed to ob- 
served Rf — J colors. This minimizes the impact of the 
poorly recalibrated Rp photometry on our sample: only 
232 objects are detected only in Rp; all other sources 
have either native SDSS colors or synthetic colors gener- 
ated from an observed Bj — J color. 

3. SPECTROSCOPIC SAMPLE 

We assembled an extensive database of follow-up spec- 
troscopy to test for contamination (e.g., background gi- 
ants, mis-identified QSOs, etc.) and bias (e.g., systemic 
misclassification of numerous M2 dwarfs as rarer M6 
dwarfs due to dim, noisy photometry) in our photomet- 
rically selected sample of low-mass stars. 

3.1. SDSS Spectroscopy 

SDSS spectroscopic observations were made with 
twin fiber- fed spectrographs, covering wavelengths from 
3800A to 9200A with a spectral resolution of A/AA ~ 
1800. Each fiber was 3" in diameter and plugged into a 
pre-drilled metal plate allowing observations across the 3 
degree field of view of the SDSS telescope. A single plate 
accomodated 640 fibers, with 320 going to each spectro- 
graph. The SDSS spectroscopic pipeline performed auto- 
mated data reduction, producing flux calibrated spectra 
corrected for telluric absorption. 

Each autumn, when the North Galactic Cap was 
completely inaccessable from Apache Point Observatory, 
SDSS-I focused its attention on an area along the celes- 
tial equator known as the Southern Equatorial Stripe, 
or Stripe 82. The observational projects conducted in 
Stripe 82 with the SDSS instruments during the autumn 
observing seasons included repeat imaging to allow in- 
vestigations of SDSS data quality and to conduct time 
variability studies, as well as stand-alone observing pro- 
grams proposed by collaboration scientists which had 
different targeting algorithms and data quality bench- 
marks than standard survey operations. Data products 



generated by these observing programs were released 
as part of DR4, and are described in more detail by 
lAdelman-McCarthv eFaLl (|2006f ). 

These additional datasets motivated us to place our 
calibration region in Stripe 82, as they contain thou- 
sands of stars with spectroscopic observations. Particu- 
larly useful to us were the SLoMaSS an d 'Spectra of Ev- 
erythi ng' sa mples, described in f ull by iBochanski et al.l 
pool and iVanden Berk et al.l (|2005l) . We identified 
stars in our sample with SDSS spectra using a 2" match- 
ing radius between the position of each spectral fiber and 
the 2MASS position of each photometric object. Due 
largely to the two observing programs mentioned above, 
10,784 stars in the calibration region, or more than a 
third of our matched sample, have an SDSS spectrum. 

3.2. Complete Observations 

For two areas within the calibration region, covering 
one square degree in total, we acquired complete spec- 
troscopic samples of low-mass stellar candidates. We re- 
fer to these areas as 'complete region 1' (4.19 < a < 
5.32, -1.04 < S < -0.58) and 'complete region 2' (11.68 
< a < 12.31, -0.24 < S < -^0.8), and the union of the 
two as the 'main complete' sample. The 'main complete' 
samples include candidate late type dwarfs satisfying the 
following criteria: 

• 2MASS astrometric positions within the bound- 
aries of complete region 1 or 2; 

• 2MASS photometry meeting the criteria expressed 
in g2XIl and J < 16; 

• SDSS or GSC counterparts meeting the selection 
criteria outlined in i i2.1.2l or I2.1.3( 

• Identification as a candidate late- type dwarf by one 
of the following two color-color cuts: 

— SDSS counterparts satisfying an r — i > 0.6 — 
1.9(i — z) color-cut, encompassing the typ- 
ical colo rs of late- t ype dwarfs identified by 
rHawlev eraII(|2002f ): 

— Objects lacking a high-quality SDSS coun- 
terpart but identified as candidate late-type 
dwarfs with a synthetic i — J > 1.2 color cut. 
These candidates were only targeted in com- 
plete region 2, where the GSC photometry un- 
derpinning the synthetic i — J color is most 
reliable. 

These cuts identified 536 targets for spectroscopic 
follow-up from an initial sample of 978 2MASS objects 
within the complete regions. The bulk of these objects 
(503) had SDSS counterparts meeting the criteria out- 
lined in i j2.1.2l for inclusion in this sample; a smaller num- 
ber (33) had photometric flaws in SDSS, and we therefore 
resorted to synthethic, GSC-based magnitudes to char- 
acterize them. 

To increase the statistics of this sample at the reddest 
colors, we expanded our criteria to cover larger areas 
of the sky, targeting candidates meeting the criteria for 
inclusion in 'extended red' and 'super red' samples. Can- 
didates in these samples met the following criteria: 
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• 2MASS photometry meeting the criteria expressed 
in g2XIl 

• Color and astrometric cuts: 

— the 'extended red' sample contains objects 
with 2MASS positions within -1.25 <S < 1.25 
and 4 < a < 6.5 or 11.5 < a < 14, and i — z > 
0.9, i- J> 2.4and J < 15.8. 

— the 'super red' sample contains objects with 
2MASS positions within the area defined by 3 
< a < 15 and -1.25 < S < 1.25, and i - J > 
2.8 and J < 16. 

The union of the 'main complete', 'extended red', and 
'super red' samples is a set of 672 candidate dwarf stars, 
of which 131 required spectroscopic observations with the 
ARC 3.5m telescope at APO. Table [1] documents the 
construction of each of these samples, listing the number 
of 2MASS sources that fell within the boundaries of each 
complete sample, the number that matched to SDSS and 
GSC sources, along with the subset that met the criteria 
for spectroscopic observation, and finally the number of 
2MASS point sources with no point source counterpart 
in either SDSS or GSC. Table [2] breaks down the sources 
of the spectra in each sub-sample. 

Finally, we note that SDSS spectra exist for many 
stars in Stripe 82 outside the boundaries of the complete 
regions, and 60 ARC 3.5m spectra were also obtained 
prior to the definition of the boundaries of the com- 
plete regions. We include these incompletely (though still 
densely) sampled sources in Table [2] under the 'Other' 
category, and used them to provide additional constraints 
on the quality of the photometric sample. 

Using the Dual Imaging Spectrograph (DIS) on the 
ARC 3.5 meter telescope at Apache Point Observatory, 
we obtained spectra of 191 candidate late type dwarfs in 
the calibration region that lacked SDSS spectra. The DIS 
spectrograph simultaneously and independently records 
spectra at blue and red optical wavelengths using a 
dichroic beamsplitter centered at 5550 A. Our observa- 
tions used the 'low' blue and 'medium' red gratings and 
a 1.5 arcsecond (3 pixel) slit, providing wavelength cov- 
erage from 3800 A to 8700 A with a typical spectral res- 
olution of ~ 700 in the blue and ~ 1000 in the red. Each 
DIS spectrum was reduced with a reduction script writ- 
ten in Pyraf, the Python-based command language for 
the Image Reduction and Analysis Facility^^. All spec- 
tra were trimmed, overscan and bias corrected, cleaned of 
cosmic rays, fiat fielded, extracted, dispersion corrected, 
and fiux calibrated using standard IRAF tasks. 

4. QUALITY CHECKS 

Measuring the luminosity and mass functions of low- 
mass stars from a purely photometric catalog requires a 
sample of the utmost quality. In particular, past efforts 
have often been waylaid by the effects of incompleteness 
(failing to detect all stars within the adopted complete- 
ness limit), contamination (incorporating objects other 

^^ PyRAF is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA. IRAF is distributed by the 
National Optical Astronomy Observatories, which are operated by 
the Association of Universities for Research in Astronomy, Inc., un- 
der cooperative agreement with the National Science Foundation. 



than main sequence late-type dwarfs in the photometri- 
cally selected sample), and bias (color-dependent errors 
in detection efficiency or source classification). We used 
our photometric and spectroscopic observations to con- 
strain the impact of these effects on our sample. 

4.1. Completeness 

Incompleteness in our sample would lead to a system- 
atic underestimate of the true physical density of stars 
in the Galaxy. Given the relative depths of the 2MASS 
and SDSS surveys, and the optical/NIR colors typical of 
low-mass stars, the inclusion of stars in our final sam- 
ple depends primarily on the 2MASS sensitivity at faint 
magnitudes. This expectation is borne out by the re- 
sults of the catalog matching described in iJ2.2l the vast 
majority ( > 99%) of 2MASS detections match to an op- 
tical counterpart. The few which do not match, however, 
represent a second possible source of incompleteness in 
our final matched sample. Below we describe tests of the 
completeness of our sample, performed by empirically de- 
riving the faint limit of 2MASS in the calibration region 
and studying the population of unmatched NIR sources. 

4.1.1. The Completeness Limit of the 2MASS survey 

Th e canonical 2MAS S 99% completeness limit is J = 
15.9 (|Cutri et al.l 120071 ) . applicable over the entirety of 
the 2MASS catalog except in areas where source con- 
fusion is the primary cause of non-detections. As the 
calibration region probes high ( |&| > 50) Galactic lat- 
itudes, stellar density is low enough that observational 
sensitivity is the main source of incompleteness. Varia- 
tions in seeing, airglow intensity, and system zeropoints, 
however, result in individual 2MASS tiles possessing true 
completeness limits up to a magnitude fainter than the 
canonical value. Deeper SDSS observations, however, 
allow a direct measurement of the magnitude at which 
2MASS becomes incomplete. 

To test the completeness limit of 2MASS in the cal- 
ibration region, added an estimate of each star's i — J 
color, derived from a least squares quadratic fit to the 
median i — J vs. i — z relation for sources with both 
SDSS and 2MASS detections, to their measured i mag- 
nitudes, producing a synthetic J magnitude {Jsyn) for 
each SDSS star. Comparing Jsy„ to the actual J magni- 
tude observed from bona fide 2MASS detections revealed 
that Jsyn is accurate to within ^^ 0.1 magnitudes (median 
error — 0.05 magnitudes, cr =0.13 magnitudes). 

Figure [3] compares the distribution of J magnitudes 
measured in the matched 2MASS/SDSS sample with 
the distribution of Jsyn values estimated for all SDSS 
sources in the calibration region, demonstrating that the 
matched 2MASS/SDSS sample becomes incomplete at 
J — 16.3. Given the ~ 0.1 mag. uncertainties in Jsym we 
select J = 16.2 as a conservative estimate of the 2MASS 
completeness limit in the calibration region. 

Applying a J <16.2 completeness limit to our cata- 
log reduced the number of 2MASS stellar candidates to 
26,585. Of these candidates, 20,869 matched to SDSS 
counterparts and 23,099 matched to GSC counterparts, 
with 18,012 candidates matched to counterparts in both 
catalouges. 

4.1.2. Understanding Unmatched 2MASS sources 
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TABLE 1 

Selecting the Complete Samples 



Sample 


2MASS sources 
in footprint 


SDSS 
Total 


counterpart | GSC counterpart | No optical 
Spectra Total Spectra counterpart 


Main Complete 
Extended Red 
Super Red 


978 
9654 
20499 


765 
7382 
15957 


503 
101 
32 


201 33 
2182 3 
4246 


12 
90 
296 



TABLE 2 
Spectroscopic Observations 



Sample 
Sample 


SLoMaSS 
Spectra 


Misc. SDSS 
Spectra 


APO 
Spectra 


Visually 
Rejected 


Main Complete 
Extended Red 
Super Red 
Other 


421 

73 

7 

7354 


15 

12 

10 

2892 


100 
19 
12 
60 




3 




TOTAL 


7855 


2929 


191 


3 



TABLE 3 

SDSS Flags For Counterparts of Unmatched 

2MASS Sources 
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Fig. 3. — Top: The number of 2MASS sources detected in the 
calibration region (solid black line) vs. J magnitude compared to 
the number predicted on the basis of SDSS detectio ns (dotted 
line). Bottom: As above, focusing on the faint limit of the 2MASS 
survey. Indicated are the adopted J =16.2 2MASS completeness 
limit (dot-dashed line), and the number of 2MASS sources matched 
to optical coun terparts (dashed line). The presence of faint 2MASS 
'orphans' is revealed by the gap between the solid and dashed lines; 
these sources are discussed in i|4.1.2l 

The 753 2MASS sources in the caUbration region with 
J< 16.3 and no optical counterpart represent a second 
potential source of incompleteness. To understand the 
nature of these unmatched sources, we visually inspected 
their SDSS and 2MASS imaging. As noted previously in 
§2.2( 550 were galaxies; these objects correctly had no 
counterpart in our optical catalogs of point sources. Of 
the remaining 203, only 11 appear to be spurious 2MASS 
detections: nine due to the mistaken identification of a 
bright star's diffraction spike as an independent photo- 
metric object, and two due to identifying random sky 



Stars Affected 



Photometric Flaw 



73 


SATURATED 


53 


PSF-FLUX-INTERP 


22 


dubious GALAXY classification 


17 


NOTCHECKED 


13 


not PRIMARY 


8 


PEAKCENTER 


3 


INTERPCENTER 


2 


undetected by PHOTO 


1 


satellite track 



2MASS diffraction spike 
2MASS false detections 



noise fluctuations as objects. 

The remaining 192 objects appear to be genuine de- 
tections of point sources, but whose SDSS counterpart 
does not meet the quality cuts discussed in §2.1.21 Table 
[3] summarizes the photometric flaws preventing 2MASS 
sources from matching to optical counterparts in the 
SDSS catalog. We exhaustively investigated these un- 
matched stars, and provide a brief summary of the most 
relevant details below. Curious reade rs can find a full 
description of this inquiry presented in ICovevI (j2006[ ) . 

SDSS counterparts for 73 of the 192 unmatched stars 
were eliminated from our sample due to the SATU- 
RATED flag being set; counter-intuitively, these were 
typically faint stars, but located within the PSF wings of 
a saturated star. The bulk of the remaining unmatched 
2MASS sources possessed SDSS counterparts with sub- 
standard imaging, either due to random noise (cosmic 
rays) or instrumental effects (landing on a bad CCD col- 
umn, etc). These flaws, indicated by the PSF-FLUX- 
INTERP, NOTCHECKED, PEAKCENTER, and IN- 
TERPCENTER flags, accounted for another 81 of the 
192 2MASS stars unable to match to their SDSS coun- 
terpart. The remaining 38 sources were excluded ei- 
ther because PHOTO identified them as extended ob- 
jects (often affecting members of a marginally resolved 
visual binary), because they were subject to extremely 
rare photometric errors (sources eluding identification by 
PHOTO, or contamination of the star's PSF by the path 
of a satellite). 

The Guide Star Catalog, an independent set of optical 
measurements, should protect against the loss of stars 
from the sample due to SDSS photometric errors. Sub- 
standard DR4 photometry indicates that the majority 
of these sources (114/186) are too dim (r < 18.5) to be 
properly detected in the Guide Star Catalog. The re- 
maining 72 sources are typically bona fide bright stars or 
close binaries with a marginally saturated component - 
the saturation contaminates both SDSS detections, while 
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the pair is unresolved in GSC imaging. 

The tendency for unmatched 2MASS sources to pos- 
sess faint and red counterparts in substandard SDSS 
photometry suggests our matching algorithm may be bi- 
ased against the latest type stars. To quantify this ef- 
fect, we compared i — J colors for sources with substan- 
dard SDSS detections to those of the full SDSS/2MASS 
matched sample. Sources with substandard SDSS detec- 
tions are skewed towards redder i — J colors than the full 
SDSS/2MASS sample, but represent less than 6% of aU 
sources for all but the reddest colors (« — J ~ 2.6, where 
small number statistics [nsta^-s/bin^S] dominate). 

From this analysis, we conclude that the sample is 
more than 99% complete to J = 16.2, and more than 94% 
complete for every color where i — J < 2.6. Photometric 
flaws in SDSS imaging marginally affect the complete- 
ness of faint sources without GSC counterparts; coupled 
with the NIR magnitude limit of this sample, this ef- 
fect may introduce a small color-dependent bias into our 
sample, such that we underestimate the true density of 
the reddest stars at the 20% level. 

4.2. Spectroscopic Quality Tests 

Using the spectroscopic catalog described in Section 
m we placed empirical limits on the contamination of 
our photometric catalog by objects other than low-mass 
dwarf stars. Additionally, we tested for bias within our 
sample, ensuring that photometric colors accurately pre- 
dict spectral types (and thus luminosities and masses) of 
late type dwarfs in our catalog. Unresolved binary sys- 
tems represent a particularly important source of bias 
for this study, and as such we discuss them separately in 

4.2.1. Identifying Exotic Contaminants 

We began our spectroscopic analysis by assigning spec- 
tral types to each spectrum, using the 'Hammer' spectral 
typing software. This set of IDL routines uses 28 spectral 
indices to estimate the spectral type of an input spec- 
trum, and then allows the user to modify the assigned 
type via visual comparison to a grid of dwarf standards. 
The automated spectral typing algorithm has been tested 
with template spectra degraded to S/N ^ 4, returning 
results accurate to ± 2 subclasses. Tests of spectral types 
interactively assigned by multiple users agree to within 
one subclass. For a full description of the Hammer al- 

prithm, we refer readers to appendix A of ICovev et al.l 

2007)12. 

Using the Hammer, we assigned spectral types to the 
vast majority of the spectra in our sample. Of the 669 
spectra in the Main Complete, Extended Red, and Su- 
per Red spectroscopic samples, only one object could not 
be confidently assigned a stellar spectral type; it con- 
tained spectral features from both a white dwarf and 
an M dwarf component. These unresolved white dwarf- 
M dwarf (WDMD) binary systems have proven to be 
relatively abundant in the SDSS spectroscopic database 
(|Silvestri et al.ll2006f ). 

^^ The Hammer has been made available for 
community use and can be downloaded from 
Ettp: //www. cf a. harvard. edu/$\sim$kcovey/thehammer , and/or 
from the tarhle provided bv iCovev et al.i (|2U07i ) to the AJ to be 
archived with the electronic edition of their article 
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Fig. 4. — Assigned spectral types as a function of photometric i — 
J color (either directly measured from SDSS/2MASS photometry, 
or transformed from GSC/2MASS photometry). Grey diamonds 
denote all calibration region objects with spectra and J < 16.2, 
while black circles indicate members of the complete spectroscopic 
samples. The dashed line demonstrates a cut to identify sources 
whose photometry and spectral types disagree significantly. 

Similarly, of the 9649 sparsely sampled spectra with 
counterparts brighter than J — 16.2, all but 24 were 
assigned stellar spectral types. Of the 24 non-stellar 
spectra, three were too noisy to confidently estimate 
a spectral type, and seven revealed additional WDMD 
pairs. The remaining 14 spectra revealed more exotic 
contaminants: four carbon stars, six unresolved galaxy 
populations (i.e, galactic bulges or clusters), and four 
QSOs. The QSOs, as weU as the seven WDMD pairs, 
fall within the u — g < 1 and g — r < 1.1 region that con- 
tains more than 95% of the QSOs in the SDSS Quasar 
Catalog (?). The galactic populations and carbon stars 
lie within the stellar locus, unidentifiable by photome- 
try alone. Extending our analysis to the 1326 additional 
spectra fainter than J = 16.2 revealed 11 more objects 
without stellar spectral types. The majority (nine) of 
these objects were QSOs, while two objects were simply 
too noisy to confidently assign a spectral type. 

In total, analysis of the 10,975 spectra within our sam- 
ple revealed 35 objects which cannot be assigned a stel- 
lar spectral type, implying a global contamination rate 
of 0.3% for our matched sample. More than half (20/35) 
of these exotic contaminants, however, can be eliminated 
by restricting our sample to objects outside the u — g < 
1 and g — r < 1.1 region. Similarly, llvezic et al.l ()2002f ) 
identified a, g — i vs. i — Ks color-color cut which can also 
help to distinguish stars from QSOs. 

4.2.2. i — J2MASS OS a Spectral Type Indicator 

Figure [4] displays spectral types of stars in our sample 
as a function of i — J color. Consistent with previous 
studies (jHawlev et al.l 120021: 1 West et al.ll2005[ ). we found 
that i — J predicts spectral type reliably in the M and K 
spectral classes. The typical spread of i — J at a given 
spectral type is ~ 0.1 mag, though there is a slight color 
dependence; the standard deviation of i — J increases 
from 0.07 mag for G5 stars to 0.12 mag for M3 stars. 

While the vast majority of the sample shows a well- 
behaved color vs. spectral type relation, there are stars 
with i — J colors significantly redder than other stars of 
the same spectral type. For our purposes, stars with late 
K and early M spectral types but i — J colors typical of 
late M and early L types are of the most concern. With 
only a small number of bona-fide late M and early L type 
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objects in our sample, consistently niisclassifying even a 
small fraction of the earlier type stars could significantly 
inflate the luminosity and mass functions at the lowest 
masses. 

To investigate the cause of these color/spectral type 
discrepancies, we inspected 20 sources with i — J colors 
significantly redder than other stars of the same spectral 
type (this cut is shown as a dashed line in Figure |4]) . 
Four of these sources were early type F and G stars that 
saturated the SDSS, and their anomalous i — J colors, 
more typical of late K/early M stars, are likely due to the 
larger errors of synthetic i — J colors calculated from GSC 
photometry. Given the large number of bona fide late- 
K/early-M objects in our sample, and the low incidence 
of mis-classified saturated SDSS stars (4/10,940 stars in 
the complete spectroscopic sample), this effect will have 
a negligible impact on our analysis. 

The remaining 16 sources with anomalous i — J 
colors are due to difficulties in properly associating 
SDSS and 2MASS sources into a matched catalog. 
One such mismatch is SMSS (Sloan M Star Survey) 
J003716. 5+000106.4; with the PSF-FLUX-INTERP flag 
set in the SDSS detection of this object, our algorithm 
incorrectly identifies a nearby faint star as this object's 
SDSS counterpart, resulting in a non-physical matched 
detection with an extremely red i — J color. Due to 
the differing spatial resolutions of the two surveys, vi- 
sual binaries with small (< 1.5") separations were of- 
ten resolved into distinct objects by SDSS, but not by 
2MASS. This scenario also results in an anamolously red 
i — J color, as the i band flux is being derived from a 
single component of the system, while the J band flux is 
the sum of both components. 

Given the relative scarcity of objects with true i — J 
colors > 3, these mismatches could artificially double the 
number of objects in the reddest bins of our luminos- 
ity and mass functions. We eliminated such mismatches 
from our sample, however, by requiring that the SDSS 
and SDSS/2MASS properties of each detection are self- 
consistent. Figure [5] shows a simple i — z vs. i — J color 
cut {i-z <i- J*0.574 - 0.738; i - 2; < 1.5) that we used 
to remove these mismatches from our sample. 

4.2.3. Kinematically Selected Subdwarf Candidates 

We estimated the number of subdwarfs in our photo- 
metric sample by identifying objects with proper motions 
indicative of halo kinematics. Such objects can be found 
photometrically using the reduced proper motion. 



Hr = r + 5 + dlogjj, 



(5) 



for a star with magnitude r and proper motion /i, ex- 
pressed in arc seconds per year. Preliminary proper mo- 
tions measured from the USNO-B and SDSS catalogs 
produce the reduced proper motion diagram shown in 
Figure [5] for objects in the calibration region with SDSS 
counterparts and > 3 a detections of proper motion 
(corresponding to a proper motion limit of 1 "/century; 
iMunn et al.ll2004) . or more than a third of the objects in 
the photometric sample (12022/30449). The disk pop- 
ulation dominates this diagram, but a spur of objects 
with reduced proper motions lying significantly below the 
disk locus is clearly seen. We isolated candidate subd- 
warfs within this spur as objects with {r — i) > 0.1 and 
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Fig. 5. — Mismatched SDSS and 2MASS detections in j-^ vs. i- 
J color-color space. Black points and greyscale contours show the 
colors of the entire photometric sample. Sources identified in the 
full spectroscopic sample as possessing discrepant i — J colors and 
spectral types are shown with nested black and white diamonds. 
Sources with discrepant i — J colors and spectral types that fall 
in the middle of the i — 2 vs. i — J stellar locus typically lack 
reliable SDSS detections, and have synthetic i — J and i — z colors 
estimated from GSC photometry and are thus consistent with the 
stellar locus by construction. Sources with i — J colors redward 
of the dashed line are likely mismatches, and are shown as filled 
circles. 




r-1 

Fig. 6. — Reduced proper motion (Hr) vs. r — i for calibration 
region stars with SDSS colors and ^ > l"/century. Objects with 
thin/thick disk kinematics are shown as black points and contours; 
subdwarf candidates are highlighted with open circles. The (r — i) 
> 0.1 and Hr > 13.86 + 6.666*{r — i) subdwarf selection criteria 
is shown as a dashed line. 

H^ > 13.86 + 6.666*(r-z), identifying 0.5% (69/12022) 
of SDSS/2MASS objects with well-measured proper mo- 
tions as candidate subdwarfs. 

The colors of these candidates butress their status as 
subdwarf candidates. These kinematically selected can- 
didate subdwarfs are ~ 0.1 mag. bluer u\ u — g than 
typical dwarfs with similar g — r colors, consistent with 
the U — B 'ultraviolet excess' previously obser ved for 
subdwarfs (see Bcssell & Wickramasingha Il979l Figure 
2). The 10 coolest (r — i > 0.55) candidate subdwarfs 
also show a redwa rd shift in q — r, consistent with the 
trend observed by IWest et alT (|2004f ) in their sample of 
spectroscopically selected SDSS candidate subdwarfs. 

We note that the 0.5% subdwarf contamination rate 
implied by this kinematically selected sample overesti- 
mates the fraction of subdwarfs present in our sample. 
The large space velocities of subdwarfs make them more 
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Fig. 7. — Measurements of the Ti05 and CaH2 indices defined by 
IReid et al.l ||1995). Black points and grayscale contours show the 
values measured from the subset of our sample with spectroscopic 
observ ations. The dashed line shows the boundary defined by Glzisi 
11997) to separate subdwarfs from solar metallicity stars. Objects 
identified as candidate subdwarfs by this cut are shown as asterisks. 

likely than slower moving disk dwarfs to fulfill a 3 tr cut 
on proper motion. Indeed, our 0.5% subdwarf contami- 
nation rate implies a local normalization of the disk and 
halo populations at the 200:1 level, a factor of two larger 
than the accepted value (|Reid fc Hawlevll2000D . 

4.2.4. Spectroscopically Identified Subdwarf Candidates 

M subdwarfs can also be identified by comparing the 
strength of metall icity sensiti ve spectral features, such 
as CaH and TiO (|Gizisl[l997l ). Unfortunately, no simi- 
larly broad features exist to identify G and K subdwarfs 
with spectra of moderate resolution and signal-to-noise 
(|Morel]| |1988^. 

We identified late-type subdwarfs in our sample by 
measuring the T i05 and CaH2 indices defined by 
iReid et all (|1995f ) for our spectroscopic sample. The re- 
sulting Ti05 vs. CaH2 diagram is sho wn in Figure U\ 
along with the polynomial fit defined by Gizii (jl997l ) to 
sel ect M type subdwarfs . Det ailed abundance analyses 
by IWoolf fc WallersteTiil (12005 ) sugg est that the subd- 
warf boundary defined bv lGizisI (|1997l ) identifies late type 
stars with [Fe/H] '^ -0.5 dex or below. 

Within the full spectroscopic sample, 41 stars were 
identified as candidate subdwarfs by their Ti05 vs. 
CaH2 ratio. Visual inspection led us to discard 22 of 
these sources with noisy or flawed spectra, leaving 19 
candidate low metallicity M dwarfs. Of 546 M stars in 
the complete spectroscopic sample, only one was identi- 
fied as a candidate late type subdwarf, representing an M 
star subdwarf contamination rate of 0.18%. Similarly, of 
the 7593 M stars in the full spectroscopic sample, the 19 
candidate late type subdwarfs represent a 0.25% contam- 
ination rate. Combining these results with the ~0.5% 
subdwarf fraction implied by the kinematic analysis in 
i i4.2.3l it appears that subdwarf/halo stars make up less 
than 0.75% percent of a matched SDSS/2MASS sample 
of late- type stars. 

4.2.5. Giant Stars 

Evolved stars, with effective temperatures and colors 
similar to main sequence stars, represent a particularly 
important source of contaminantion to a photometric 
catalog of main sequence stars. Since giants are much 
more luminous than dwarfs, a given flux limit for ob- 
servational detection implies a larger distance limit for 



giants than for main sequence dwarfs of the same color; 
effective distance limits for giants are typically 50 times 
larger, such that a survey is sensitive to giants over a 
volume more than 100,000 times greater than that for 
dwarf stars. Thus, while the absolute space density of 
giants is low, they can make up a significant component 
o f a magnitude lim i ted sa mple o f point sources. 

iBessell fc Bret^ (Il98l and iWorthev fc Le3 (l2006ll 
have identified the JHKg and VIKg color spaces (see 
their Figs. 5 and 4, respectively) as particularly use- 
ful for separating late- type dwarfs and giants. Given 
the size of the expected JHK color shifts (~0.2 mag) 
relative to the accuracy of the 2MASS photometry, we 
can only draw the general conclusion that the sample is 
dominated by dwarf stars; in our sample, J — H becomes 
noticably bluer redward oi H — Kg ~ 0.15, with no dis- 
tinct plume of giant stars extending towards J — H ^^ 
1 a.t H — Kg ~ 0.15. Very few stars in the calibration 
region have giKg colors indicative of giants, and the tens 
of objects that do make up a very small percentage of the 
total number of sources detected at that i — Kg color. 

Spectroscopic indices allowed us to identify G and 
K giant candidates within our sample, producing an 
empirical estimate of their effect on our analysis. 
Using 2388 spectra from muncrous spectral libraries 
(Torr es- Dodgen fc W e aver 1 993; Danks fc D en nefeU 
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tested spectroscopic indices that identify giant stars in 
moderate resol ution, moderate signal to noise spectra 
(jMalvuto et al.l 119971 : iMalvuto fc Schmidt-KaTe^ |l999). 
These indices are given in Table |4j with each index 
calculated as a ratio of the mean flux in the spectral 
regions identified as the index numerator and denomi- 
nator. The only exception is the Mg2 index, which was 
calculated as a psuedo- equivalent width , consi stent with 
its original definition in iMorrison et all (|2003f ). 

We used Mg2 vs. g — r as our primary giant/dwarf 
discriminant (Morrison et al., in prep; Yanny et al., in 
prep.). At red g — r colors, however, the dwarf and giant 
sequences reconnect. To identify even the reddest giant 
stars, we selected giant candidates from this confused re- 
gion of Mg2 vs. g — r space using three supplementary 
spectroscopic criteria (Na D 5900 vs. Mg b; TiO 6700 
vs. Blend 6497; CN 7900 vs. Na D 8200) to reject likely 
late type dwarfs. This technique identified 104 candi- 
date giants from 8750 objects in our sample with both 
g — r colors and spectra, representing a 1% global giant 
contamination rate. 

To test if these spectroscopically identified candidate 
giants are consistent with standard models of the stel- 
lar population of the Milky Way, we generated simulated 
SDSS/2MASS obser vations using the TRILEGAL code 
(jGirardiet al.ll2005f ). The top two panels of Figure [8] 
show the locations of giant stars in i — J vs. J color- 
magnitude space as predicted by the TRILEGAL sim- 
ulation using standard Galactic parameters and the lo- 
cation of the calibration region (top left panel) and as 
identified in our spectroscopic sample (top right panel). 
The TRILEGAL simulation predicts that giants should 
reside in a relatively narrow locus in i — J vs. J color- 
magnitude space stretching from J ~ 6 and i — J ^ 
1.3 down to J ~ 16 and i — J ^ 1.0. The locations 
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TABLE 4 
Spectral indices for Dwarf/Giant discrimination 



Band name 


Numerator 1 


Numerator 2 


Denominator 


Reference 


Mg2 


4935- 


-4975 


5303- 


-5367 


5130 - 


-5200 


Morrison et al. (20031 


Mgb 


5100- 


-5150 


5320- 


-5420 


5150- 


-5320 


Malvuto et al. (19971 


Na D 5900 


5789- 


-5839 


6020- 


-6120 


5839- 


-6020 


Malvuto et al. (1997) 


Blend 6497 


6406- 


-6470 


6582- 


-6637 


6470- 


-6536 


Malyuto & Schmidt-Kalcr (19991 


TiO 6700 


6425- 


-6525 


6970- 


-7040 


6613- 


-6832 


Malvuto ct al. (19971 


CN 7900 


7815- 


-7887 


7970- 


-8037 


7887- 


-7970 


Malvuto & Schmidt-Kaler (1999) 


Na D 8200 


8120- 


-8170 


8210- 


-8260 


8170- 


-8210 





of our spectroscopically identified candidate giant stars 
in i — J vs. J color-magnitude space match tliese pre- 
dictions well; though the spectroscopic catalog only con- 
tains sources with relatively faint J magnitudes, Figure 
[5] shows the faint spectroscopically identified giant star 
candidates are indeed concentrated about i — J '^ 1.1. 
Additionally, the photometric sample shows a distinct 
plume of sources tracking the expected giant star locus, 
most prominently at J ~ 11 and i — J ~ 1.2. 

The bottom panels of Figure [5] show histograms of the 
i — J colors of all stars, as well as just giants, in both 
the TRILEGAL simulation and the observed data. The 
TRILEGAL simulation predicts a global contamination 
rate for our sample of 4%, but this is a strong func- 
tion of color, with giant star contamination rising above 
10% blueward of i — J ~ 1.15. We drastically reduced 
giant star contamination by restricting our analysis to 
sources with J > 12 and i — J > 1.25 (> K5), bounded 
in Figure[8]by blue dashed lines. These cuts reduce giant 
star contamination to a negligible level (less than 2% in 
the bluest, most contaminated color bin) while sacrificing 
only a handful of bona fide low-mass dwarf stars. 

5. THE LUMINOSITY FUNCTION 

The previous section demonstrates that we can con- 
struct a matched SDSS/2MASS dataset for stars with 
12 < J <16.2 and 1.25 < i — J < 4.5 where incomplete- 
ness, contamination, and bias are limited to <1% for the 
full sample, and <5% for stars in a small (~ 0.1 mag) 
color range. We are thus confident in the quality of the 
calibration region data, and now describe the use of this 
dataset to measure the luminosity and mass functions of 
the Galactic disk. 

Begining with the combined photometric sample as- 
sembled in 321 we culled the catalog to ensure complete- 
ness and reduce contamination and bias from background 
giants and bad matches. In doing so, we restricted our 
analysis to only those sources meeting the following ad- 
ditional color and magnitude criteria: 

• completeness cut: J < 16.2 & SDSS or GSC optical 
counterpart (leaves 25,947 candidates); 

• giant contamination cut: J > 12&z — J> 1.25 
(leaves 13,159 candidates); 

• bad match cut: i — z < 1.5 & i — z > 0.574(i — J) 
- 0.738 (leaves 13,088 candidates) 



QSO/WDMD contamination cuts 



13. 



5 > 1 or 



g ~ r > 1.1 or i — K < 1.5 -|- 0.8 *(<? — «) (leaves 
13,064 objects in final sample). 

We will refer to this set of 13,064 objects as the Lu- 
minosity Function (or LF) sample. All but 714 of these 
objects possess native SDSS magnitudes; the rest possess 
psuedo-SDSS magnitudes transformed from GSC obser- 
vations. 

The calibration region probes a sight line that passes 
directly through the mid-plane of the Milky Way, where 
the bulk of the dust in the Galaxy is thought to lie 
(lidust ^ 125 pc; Marshall et al. 2006). As the vast 
majority of the stars in our sample lie at distances in 
excess of 250 pes, or two dust scale heights past the mid- 
plane of the Galaxy, we assume the bulk of the extinc- 
tion along the line of sight to each star is foreground 
extinction. We therefore applied extinction corrections 
to the photometry for stars in our sample, using the ex- 
tinction estimates measured byiSchlcgel ct al. (1998) and 
stored in the SDSS database (for GSC detections, we 
adopt the reddening estimate of the nearest SDSS star 
in our sample). We use the relative extinctions in Ta- 
ble 6 of ISchle gel et al.l ()1998l ) to apply corrections to our 
SDSS and 2 MASS photoni e try us ing t he Rv=3.1 extinc- 
tion laws of iCardelli et~aI1 (pMl and lO'DonneTll (fT99l . 
These corrections are minimal, with a maximum E(z — J) 
of 0.11 mags and more than 95% of the sample possessing 
E(i — J) < 0.05 mags; indeed, tests indicate that neglect- 
ing the effects of extinction produces a negligible impact 
on our results. 

5.1. The SDSS/2MASS Color- Magnitude Relation For 
Late-Type Stars 

Producing an LF from our photometric, magnitude- 
limited sample required us to estimate the absolute mag- 
nitude and distance of each star from its observed colors 
and an adopted color-magnitude relation (CMR). The 
CMR for low-mass stars is reasonably well known in the 
standard Johnsons-Cousins photometric system; the rel- 
ative youth of the SDSS, however, make the transfor- 
mation of standard CMRs onto the SDSS photometric 
system difficult for stars with red colors. 

By spectroscopically classifying a sample of low- 
mass st ars wi th me asured SDSS and 2MASS colors, 
iHawlev et al.l (|2002[ ) derived photometric parallax re- 
lations for low-mass stars in SDSS by linking distinct 
color-spectral type and absolute r nagnitude-spectral 
Subsequent studies (West et al.l l2005t 



type relations. 



^^ The need for u — g colors to implement these cuts restricts 
their usage to only those objects with SDSS optical counterparts. 



iBochanski et al.ll2007tlWest et al.ll2008l1 have used larger 
samples to produce new measurements of the relation be- 
tween spectral type and SDSS color, shifting the mean 
colors of a given spectral type somewhat bluer; these new 



12 



Covey et al. 



A I I I I I I I I I I I I I I "1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



10 



12 



14 



a 



16 

1000 

100 



10 




Fig. 8. — The location of giant star candidates in j — J vs. J color-magnitude space, as predicted by the standard TRILEGAL galaxy 
model (left panels) and as indicated by our photometric and spectroscopic catalogs (right panels). Dwarfs are shown as black points and 
contours in the top two panels, while giant stars are highlighted as red asterisks and identified using surface gravity in the TRILEGAL 
output or spectroscopic indices in the SDSS/2MASS sample. Dashed blue lines show photometric cuts of J > 12 and i — J >1.25 (> K5). 
Bottom panels comp are the i — J color distributions for all stars with J > 12 (solid black line) and for giants (dashed red lines). To reduce 
contamination by giants, we restrict our luminosity function sample to sources with J > 12 and i — J > 1.25 (blue dashed lines). 



color-spectral type relations produce new photometric 
parallax relations when coupled to the previously mea- 
sured absolute magnitude-spectral type relation. Stellar 
colors and magnitudes in the SDSS and 2MASS systems 
have also been predicted usin g theoretical s t ellar mod- 
els, such as t hose calculated bv lGirardi et al.l (|2002D and 
iDotter et all (|2007il . 

Ultimately, it would be best to determine the CMR 
by directly measuring native SDSS colors and magni- 
tudes for stars with measured trigonometric parallaxes. 
In practice, nearly all such stars saturate the SDSS 2.5m 
camera, requiring observations with smaller telescopes 
to be transformed o nto t he 2.5m photome t ric sy stem. 
IWilliams etHI (|2002[ ) and lGolimowski et all (|2008f ) have 
observed parallax standards with the 0.6 meter SDSS 
Photometric Telescope (PT) and the USNO 1 met er tele- 
scope, using the relevant tra nsformations (Tucke r et al.l 
120061 : iDavenport et al.l l2007t ) to place their photometry 



on the SDSS 2.5 meter system. We fit preliminary pho- 
tometric parallax relations to describe the dependence 
of Mj on r - J and i - J using 2MASS JHK mag- 
nitudes and tran sformed SDSS uqr iz magnitudes kindly 
provided to us bv lGolimowski et al.l ([200S) prior to publi- 
cation^'*. These empirical CMR fits are shown in Figure 
[9l along with SDSS/2MASS colors and magnitudes for 
individual stars measured in Appendix A. 

The empirical CMR we have derived, however, suffers 
from some level of uncertainty. In particular, the data 
are too sparsely sampled to definitively reveal the pres- 
ence or absence of a sharp break at i — J ^ 2.3, as is 
seen in other color-magnitude diag rams on the J o hnson - 
Cousins system. As discussed by iReid fc GizisI (|1997D , 

^■^ With an ultimate goal to measure the MF of low-mass 
stars, we chose to use CMRs to estimate Mj, for which well- 
established mass/luminos ity relations exist IHenrv & McCarthy! 
[T993: [Delfosse et al.|[2000l) . 
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Fig. 9. — Color-magnitude relations necessary to estimate the 
distribution of stellar absolute magnitudes from the i — J colors 
of low-mass stars in the calibr atio n region. The solid line shows a 
preliminary fit to the data of IGolimowski ct al. (2008), assuming 
a smooth fit through j — J ~ 2.3; the dotted line shows a fit to 
the same data, but assuming the presence of a sharp break at 
i — J ~ 2.375. Shown for comparision (but not fitted directly) 
are the SDSS/2MASS colors and magnitudes of individual stars as 
reported in Appendix A. 

this feature significantly affects the shape of the LF in- 
ferred from a magnitude limited photometric sample. In 
the absence of data clearly indicating the presence of a 
sharp break in the SDSS/2MASS CMR, for our analysis 
we adopt the smoothly varying CMR shown as a solid 
line in Fig. [9l Understanding the important role the 
CMR plays in measuring the luminosity and mass func- 
tions of magnitude limited samples, however, we include 
uncertainty in the empirical CMR as a potential source 
of systematic error in our analysis, which we discuss in 
more detail in i)5.4l 



5.2. The Veff Technique 

Robust space density measurements of low-mass stars 
required a careful assessment of the volume we sampled 
for stars at each luminosity. In particular, we needed 
to account for the interplay between the vertical density 
gradient of the Galactic disk and the magnitude limits 
of our catalog. We did so by using the effective volume 
technique. 

With roots in the y^ technique first introduced by 
ISchmidtl (|1968D . t he concept of an 'effective volume' was 
first suggested bv iHuchra fc SargentI ()1973D . who modi- 
fied the y^ test to accomodate a spatial gradient as a 
function of position on the sky. Felted (jl976f ) extended 
the Ve// formalism to include gradients as a function of 
depth along the line of sight and proved that the sum of 
y^ for a sample of objects can serve as an unbiased es- 
timator of their spatial density, even in situations where 
the objects are distributed inhomogeneously. 

The effective volume of an object can be expressed 
mathematically as 



K 



eff 



n 



dmax 



dmin 



Po 



(6) 



where fl is the observed solid angle, r is the position along 
the line of sight to an object, dmin is the minimum dis- 
tance at which the object could be detected, dmax is 
the maximum distance at which the object could be de- 
tected, and y- describes the ratio of the sampled density 
to the local density as a function of position along the 
line of sight. For non-uniform source distributions, — 
effectively weights each volume element along the line 
of sight according to the expected likelihood of finding 
a source in that volume, such that sampling equivalent 
effective volumes should result in samples with similar 
numbers of sources. 

Using this methodology, space densities (represented 
by <&) are given by the sum of p^ for each source within 
the sample: 



*=i:i7^=E 
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(7) 
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As shown bv lHuchra fc SargentI ()1973( ) , the uncertainty 
in this estimator of <& is: 



(7$ 




(8) 



Combining the bright and faint limits of our survey 
sample (J = 12 and J = 16.2 respectively) with the value 
of Mj estimated for a single source allows us to calculate 
its individual distance limits, dmin and dmax. Thus, all 
that remains in the calculation of Ve//, and thus ^{Mj), 
is to integrate r^ — between dmin and dmax. 



po 



5.3. 



Adopted Galactic Structure Model 

Adopting a standard multi-component Galactic model 
allows us to calculate the stellar density profile expected 
along a given line of sight; the total density, ptot, is sim- 
ply the sum of the density profiles of the individual disk 
and halo components. Expressed as a function of po- 
sition within the Galaxy (using the standard Galactic 
coordinates R and Z), these profiles are: 



Pthin = exp 
Pthick = exp 

Phalo 



( L^^ ( \3L.\ 

( \.IL-] ( IJ1L-] 

V Hfh.akJ gx„V ^th^ak J 



a. 



lR^ + {Z X -)2 



Ptot — /l 



thi' 
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I 
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The total stellar density at a given point in the Galaxy 
thus depends on the nine parameters given in Table [5l 
Given the Galactic latitude (& ~ -62) and the color and 
magnitude limits of the calibration region sample, — is 
set primarily by the vertical profile of the thin disk; the 
thick disk only begins to dominate observed star counts 
at distances of ~1.25 kpc, where even our most intrinsi- 
cally luminous stars approach the faint limit of our sam- 
ple. 

In principle, the parameters in Table [5] can be con- 
strained by fitting observed star counts to the output 
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TABLE 5 
Galactic Structure Parameters 



Parameter name Parameter description 



Adopted value 



fthin 


thin disk stellar density { 
in the solar neighborhood 


1 - fthick - fhalo) 


Hthin 


thin disk scale height 


280 pc 


Lthin 


thin disk scale length 


2500 pc 


fthick 


thick disk stellar density 
in the solar neighborhood 


5% 


Hthick 


thick disk scale height 


900 pc 


Lthick 


thick disk scale length 


3500 pc 


c/a 


halo flattening parameter 
axial ratio c/a 


0.7 


rhalo 


exponent of halo power 
law density gradient 


-2.75 


fhalo 


halo stellar density 

in the solar neighborhood 


0.15% 



of Galactic stellar population models produced with a 
range of Galactic structure parameters. The best es- 
timates, however, are produced by comparing models 
to star counts along multiple lines of sight that sample 
regions of the Galaxy dominated by different Galactic 
components. Even in this best case, there often remain 
significant degeneracies between various Galactic struc- 
ture parameters (e.g., fthick and Hthick; see Fig. 1 of 
iSiegel et al.l[200l) . 

As the calibration region probes only a single sight line 
through the Galaxy, attempts to fit the Galactic model 
directly from the observed star counts produced unphys- 
ical results with significant degeneracies between param- 
eters. A follow-up paper in this series (Bochanski et al. 
2009, in prep.) will present robust Galactic structure pa- 
rameters measured from low-mass stars detected across 
the entirety of the SDSS footprint; for now, however, we 
calculated — using standard Galactic structure parame- 

ters from the lite rature (see Table 1 of ISiegel et al.ll2002l : 
iJuric et"ar ''2008') which are shown in Table [51 As Figure 
[TUl demonstrates, the profile of stellar counts predicted 
by this model is in good agreement with the observed 
star counts in the calibration region over the entire mag- 
nitude range of our sample. 

5.4. The Luminosity Function of the Galactic Disk 

Applying the Ve// technique to the LF sample, we ob- 
tained a first measure of the LF of the Galactic disk, 
shown in the top panel of Fig. [11] as green filled di- 
amonds. Derived from a magnitude-limited sample of 
stars with distances up to a kiloparsec or more, this mea- 
surement suffers from Malmquist-type effects and rep- 
resents the LF of systems, not individual stars, in the 
Galactic disk. 

To account for the effects of unresolved binaries and 
Malmquist-like biases, we constructed an algorithm to 
reveal the intrinsic stellar LF required to reproduce the 
raw system LF measured from the calibration region sam- 
ple. This routine uses a Monte Carlo approach to pro- 
duce a simulated photometric sample given an input LF, 
CMR, and parameters to describe the adopted Galac- 
tic model and binary population. We then conduct a 
Ve// analysis on this synthetic sample, which incorpo- 
rates both observational errors and the intrinsic width 
of the CMR. Using the difference between the synthetic 
'measured' LFs and the actual raw LF measured from 
the calibration region to adjust the input LF used to run 
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Fig. 10. — Star counts as a function of distance for stars with 
similar Mj values (colored circles). For clarity, each observed pro- 
file is shown twice: once compared to a common star count profile 
predicted for the Galactic structure paramaters adopted here (solid 
black line, bottom), and once compared to the same model profile 
but with a unique vertical offset (dashed grey line). Observed star 
count profiles have been normalized such that the total number of 
stars along the line of sight matches the prediction of the model, 
allowing sets of stars with different local density normalizations 
to be compared to the same model profile. The Galactic struc- 
ture model adopted here produces a good agreement between the 
observed and predicted star count profiles. 

the next simulation, we iteratively converged on the true 
LF required to replicate our observed LF in the presence 
of these observational biases. 

For our primary analysis, we adopted the fiducial 
Galactic model outlined in Table \5\ We began by gener- 
ating a synthetic set of star counts, randomly sampling 
the raw calibration region LF as an initial guess for the 
input LF and accounting for the increase in Ve// towards 
brighter Mj. We used the smooth empirical CMR shown 
as a solid line in Fig. [9] and the Mj vs. mass relation 
derived bv lDelfosse et al.l (|2000f )^^ to assign masses and 
i — J colors to each star. We included an intrinsic disper- 
sion in the CMR by sampling color offsets from a normal 
distribution with cf color '^^ 0.1 mag, or equivalently, mag- 
nitude offsets with amag ^ 0.5 mag. 

We accounted for the presence of unresolved binaries 
in our sample by randomly assigning secondary com- 
panions to synthetic stars with a mass dependent bi- 
nary fraction, fbin{M) = 0.45 — (0.7 — Mp)/4, consis- 
tent with observa tions of the stellar binary fraction in 
the Galactic disk iD ugucnnoy et al. 1991; Rcid fc GizisI 
19971: iDelfosse et al.. .2004: .Lada . .2006.: .Burgasser et al.l 



20071 ). In our model, companion star luminosities were 
randomly selected from the same input LF used to gen- 
erate the primary stars, but secondaries were required 
to have lower luminosities than their primaries. This 
quasi-independent sampling of the LF nicely reproduces 
the observed correlation between mass ratio and primary 
mass, producing a mass ratio distribution which peaks at 
0.4 for stars more massive than 0.6 Mq and at 0.9 for 



15 We note that the IDelfosse et al.l I I2000I ) Mj vs. mass re- 
lations as published and archived on the A&A website and the 
arXiv preprint server contain typographical errors. The orig- 
inal authors confirm that the correct relation is Log -rj— = 

Mq 

(10-3 X (1.6 -I- e.OlMj -I- 14.888M5 - 5.3557M5 -I- 0.285181MJ). 
(T. Forveille, priv. communication, 2008) 
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Fig. 11. — The calibration region luminosity function. Top 
Panel: System luminosity function measured from the calibration 
region (green triangles), and the resultant single-star luminosity 
function, after accounting for the impact of unresolved binaries 
and Malmquist type effects (black circles). Middle Panel: Sys- 
tematic uncertainties in the luminosity function, shown by com- 
paring the calibration region single-star luminosity function (black 
circles) with results obtained when adopting a color-magnitude re- 
lation with a sharp break at ~M4 (i — J =2.375; red dashed line), a 
Galactic model with a ~20 scale height (green dash-dot line), and 
correcting for binaries assuming a fiat mass ratio distribution (blue 
dotted line). Bottom Panel: The calibration region single-star lu- 
minosity function (black circles), including error bars for statistical 
(hatted) and systematic (un-hatted) effects. Shown for compari- 
son are luminosity functions d erived from the v olume complete 8 
pc sample (blue filled squares; IReid et al.ll2002|) and for la te-M/L 
dwarfs within 20 pc (red filled triangles: ICruz et al.ll2007l) . 

stars less massive than 0.2 M©. 

We assumed that all binary systems in our sample are 
unresolved, based on the relatively poor spatial sampling 
of the 2MASS detectors (2"pixels). This assumption is 
not uniformly correct, as resolving a binary depends on 
the ratio of this resolution limit to the binary's apparent 
separation, itself a function of both the physical sepa- 
ration of the companions and the distance to the sys- 
tem. Systems with an Mj = 8 primary, for example, 
are resolved only if their separation exceeds 30-200 AU, 
depending on the location of the system within the vol- 
ume to which we are sensitive to such stars; the situation 
is even worse for instrinsically brighter stars, which are 
drawn from a more distant volume and thus have smaller 
apparent separations for equivalent physical separations. 
As most binary stars have separations less than 100 AU, 
and our sample is dominated by stars at the largest dis- 
tances, we treat all simulated binary systems as unre- 
solved. Fluxes from each companion are merged to pro- 



duce system magnitudes and colors, at which point ran- 
dom observational errors are added to produce a psuedo- 
observed photometric sample. 

We then measured the LF of this sample using the same 
algorithm as applied to our observed sample, and used 
the difference between the two measured luminosity func- 
tions to modify the input LF used to produce the syn- 
thetic sample. After a small (n~5) number of iterations, 
the LF measured from the synthetic sample converges to 
that measured from the calibration region dataset. The 
final input LF whose 'measured' LF reproduces the raw 
calibration region LF is therefore a measure of the true 
LF of the calibration region after correcting for unre- 
solved binaries and Malmquist- type biases, as each effect 
is included in producing the catalog of synthetic stars. 

The calibration region LF with corrections for 
Malmquist biases and unresolved binaries is shown with 
black circles in each panel of Fig. [11] Somewhat counter- 
intuitively, aside from making the peak at Mj — 7 more 
prominent, the corrected LF is very similar to the uncor- 
rected LF we measure from the calibration region. Ap- 
plying binary and Malmquist corrections independently 
showed that the two effects offset o ne another to a large 
degree; as discussed extensively by IStobie et alj ([1989), 
the Malmquist correction shifts the peak of the LF to 
brighter magnitudes, primarily by decreasing the den- 
sities measured on the faint side of the LF peak. Ac- 
counting for unresolved binary systems, however, has the 
opposite effect - splitting binary systems into individual 
components increases the source density in the stellar LF 
above that of the system LF, and is most important for 
fainter stars, which are most likely to be hidden as secon- 
daries to more luminous primaries. The primary impact 
of these corrections, therefore, is to increase the peak 
density of the LF (at Mj = 7.5) by 20%; outside the 7 
< Mj < 8 range, the corrected stellar LF is remarkably 
similar to the measured system LF. 

To test how sensitive our LF measurement is to our un- 
derlying assumptions, we re-derived the LF after modi- 
fying the adopted CMR, Galactic model, and binarity 
prescription. We varied, in turn, the binary prescription 
to assign secondaries assuming a flat mass ratio distribu- 
tion, the Galactic model by increasing the scale height of 
the thin disk by 50 parsecs, and the CMR by adopting the 
red dashed empirical fit in Fig. [Sjthat includes a sharp 
jump in magnitude at i — J = 2.375. The LFs derived 
under each of these modified assumptions are shown in 
the middle panel of Fig. [11] 

We characterized the systematic uncertainties in our 
analysis as the fractional change induced in the LF by 
modifying each portion of the LF measurement algo- 
rithm; adding each of these individual uncertainties in 
quadrature produced a global estimate of the system- 
atic uncertainties in our analysis, which is shown as an 
additional set of error bars for the corrected single-star 
LF in the bottom panel of Fig. [TT] The binary correc- 
tion introduces a relatively small ^^5% systematic uncer- 
tainty, rising above 10% only at the faintest magnitudes. 
The LF is somewhat more sensitive to the Galaxy model 
adopted in the Vg// analysis, typically producing a 20% 
uncertainty in the derived LF. The uncertainty in the LF 
due to the Galaxy model is more pronounced at brighter 
magnitudes, where the volume sampled is several scale 
heights above the plane of the Milky Way, and thus the 
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TABLE 6 

The Luminosity Function of the 

Calibration Region 





Raw 


Corrected 






Mj 


LF 


LF 


o"stat 


O-sys 


5.25 


2.79 


2.99 


0.25 


0.79 


5.75 


3.41 


3.82 


0.31 


0.83 


6.25 


4.27 


4.83 


0.39 


1.20 


6.75 


5.70 


6.53 


0.54 


1.30 


7.25 


8.09 


9.64 


0.77 


1.55 


7.75 


9.52 


11.33 


1.04 


1.84 


8.25 


8.40 


8.08 


1.23 


0.72 


8.75 


7.01 


6.51 


1.46 


2.82 


9.25 


5.37 


4.71 


1.67 


4.54 


9.75 


6.12 


4.60 


2.37 


1.31 


10.25 


3.26 


3.37 


2.35 


0.79 


10.75 


1.42 


1.10 


1.75 


1.27 


11.25 


0.30 


0.41 


0.59 


0.28 


Note. 


— LFs 


and uncertainties re 


ported 


in units of 10^^ 


stars (0.5Mj 


)-v- 





densities expected by the Vg// calculation are more sen- 
sitive to the adopted thin disk scale height. 

The uncertainty in the exact behavior of the CMR at 
i — J ^ 2.375, however, introduces the largest system- 
atic uncertainty into our analysis; LFs produced using 
the two empirical fits in Fig. [9] differ by more than 
90% at Mj — 9.25. Reducing this uncertainty is clearly 
necessary to realize the full promise of the SDSS and 
2MASS to characterize the low-mass stellar content of 
the Galactic disk. Doing so will require a concerted ef- 
fort to transform photometry of trigonometric parallax 
standards from moderate aperature telescopes onto the 
native photometric system of the SDSS 2.5m. 

The bottom panel of Figure [TT] displays the single-star 
calibration region LF along with error bars accounting 
for both statistical and systematic effects. Shown for 
comparison are t he lo w-mass stell a r LFs measured by 
iReid et al.l (|200l and iCruz et al.l (|2007l) from volume 
complete samples in the solar neighborhood; with the 
exception of the Mj = 6.75 bin, the calibration region LF 
agrees with the LFs measured from the volume complete 
samples, given the uncertainties in each measurement. 

6. THE MASS FUNCTION 

The same simulations that allowed us to correct our 
observed LF for unresolved binaries and Malmquist ef- 
fects also provided an opportunity to measure the stellar 
mass function. Applying the Veff technique to the syn- 
thetic stellar sample, but binning by mass instead of Mj 
magnitude, produced the MF shown in the top panel of 
Figure fT2l again with both systematic and statistical er- 
ror bars. The logarithmically binned mass function of 
the calibration region has a clear peak at M* ~0.3 M©. 
We note that the systematic error bars, due mainly to 
the uncertainty in the color magnitude relation, domi- 
nate the statistical error bars for all mass bins; until the 
SDSS/2MASS color-magnitude relation is clarified, we 
will be unable to reap the full benefit of a sample of this 
size. 

We used a Markwardt minimization to fit various an- 
alytic forms to this data, with results shown in the top 
panel of Fig. [121 We fit a log normal form to the data, 
finding best fit values of A=0.1 (in units of stars (Log 
Mq)~^ pc~^ for comparison with the values reported by 
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The 



e-star (not system) mass function (- 



^ dlogM ' 

of the Galactic disk as measured from the cahbration region sam- 
ple (aste risks; hatted error bars show statistical errors, unhatted 
error bars show syste matic uncertainties). Top Panel: Ana- 
lytic fits to the calibration reg ion MF data; log normal, power 
law, and broken power law fits to the data are sh own as dotted, 
dot-dashed, and dashed lines respectively. An F-test indicates t 
hat the log normal fit provides the greatest fidelity given the num- 
ber of parame ters used to produce the fit. Middle Panel: A 
co mparison between the calibration region MF and that derived 
by IReid et al.l II2002I) from the 8-pc volume complete sample (red 
diamonds); the two measurements agree within the uncertaint ies 
for all but the lowest mass point. Bottom Panel: Analytic descri 
ptions o f the single-star MFs presented by Reid et al. ( 2002) and 
IKroupal (2002), shown as red dotted and blue dashed lines respec- 
tively, compared to the log normal fit adopted here (solid black 
line). While the data underlying these fits agree well, the analytic 
descriptions diverge significantly, particul arly at the low-mass end. 



IChabriei|[2003[) . Mc=0.29, and cr=0.32; we also fit a sin- 
gle (a =1.1) and broken power law (Q;/ogM>-o.5=2.04; 
a;ogM<-o.5=0.2) to the data as well. An F-test evalu- 
ates the quality of an analytic fit to a dataset relative to 
the number of degrees of freedom in the fit. Applying 
an F-test to these fits revealed that the log normal for- 
malism produces a quantifiably better fit than a single 
power law: there is a less than 1% chance that the un- 
derlying distribution is truly a single power-law and that 
the improved quality of the log-normal fit is a random 
occurance. A dual power law fit provides an even better 
fit than either the single power law or the log normal 
form, but at the cost of additional degrees of freedom; 
the F-test indicates there is a 35% chance that the dual 
power law fit does not provide a statistically superior fit 
to the data, and as such we adopted the log normal fit 
as the most efficient analytic description of our data. 
Mc represents the peak of the logarithmically binned 
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TABLE 7 




The Mass Function of 


THE 


Calibration REGIO^ 


1 


Log M* 


MF p 


Cstat. 


0"sys. 


-0.15 


5.42 


0.44 


1.29 


-0.25 


6.13 


0.50 


1.53 


-0.35 


7.91 


0.65 


1.52 


-0.45 


10.76 


0.89 


1.73 


-0.55 


10.82 


1.18 


1.54 


-0.65 


7.40 


1.28 


1.37 


-0.75 


6.07 


1.59 


3.54 


-0.85 


4.73 


1.93 


3.50 


-0.95 


5.61 


3.38 


1.47 



Note. — MF densities and uncer- 
tainties are reported in units of 10^ 
stars (O.llogM. )~^pc"^. 

mass function, such that while the density of stars per 
logarithmic mass bin declines for masses less than Mc, 
the width of each bin in linear mass units is also changing 
as well. As a result, a log normal MF with a given Mc 
corresponds to a linearly binned MF whose peak occurs 
at smaller masses than Mc. The log normal fit described 
above corresponds to a linearly binned MF that turns 
over at 0.17 Mq. 

The MF we measure here agrees well with previous in- 
vestigations. Using a Monte Carlo routine to refit the 
MF after applying random offsets consistent with each 
datapoint's total error budget, we derive 90% confidence 
intervals for the log-normal fits' parameters of A=0.08- 
0.12, Mc=0.20-0.50, and q- ^^O. 22-^0.37. Th ese agree well 
with the values reported bv lChabrieil (|2005f ). with the ex- 
ception of (T, for which Chabrier finds a somewhat larger 
value of 0.55, outside the bounds of our 90% confidence 
interval. This may result from the difference in the range 
of masses studied here and by Chabrier; Chabrier's work 
attempts to describe the MF well above and below the 
masses sampled here, which may drive the fit to larger 
values of a, and wider MF peaks. We also include in 
the middle panel of Fig. [T2lthe indi v idual datapoints for 
the MF constructed by R eid et al.l (|2002f ). which agree 
within 1 a with the datapoints reported here for all but 
the lowest mass bin. Similarly, the a = 1.1 power law fit 
to the entire mass range sampled he re is a good matc h 
to the a— 1.2 power law measured bv lReid et al.l ()2002l ). 

Simply comparing analytic fits from other studies to 
the MF data measured here, however, can give the po- 
tentially false impression of disagreement in the underly- 
ing MFs. This is due in part to the sensitivity of analytic 
fits to the exact formalism chosen to represent the data, 
as well the exact range of masses spanned by the fit. As 
shown in the bottom panel of Fig. [121 the log normal de- 
scription of the MF data reported here diverges stro ngly 
from the analy tic fits reported bv lReid et al.l (J2OO20 and 
lKroupal(|2002l ). which would seem to provide evidence for 
significant MF differences. In particular, the log normal 
fit presented here predicts an MF that peaks in the stellar 
regime, while the break point s adop ted fo r the power-law 
fits measured by iReid et al.l (|2002l ) and iKroupal (|2002l ) 
imply that the MF rises all the way to the sub-stellar 
limit. As shown previously, ho wever, the individu al MF 
datapoints derived here and bv iReid et al.l (|2002l ) agree 
well, revealing that these varations are due more to dif- 
ferences in the formalisms adopted to describe the MF 



than to confident detections of significant MF differences. 

Traditionally, the utility of analytic descriptions of the 
data has been to produce a more robust description of 
a potentially noisy dataset, collapsing the many degrees 
of freedom of the sample to improve the quality of the 
measurement of a reduced set of parameters. Given the 
richness of modern-day astronomical datasets, however, 
where one can construct samples of stellar populations 
that contain thousands, and in some cases millions, of 
objects useful for deriving mass functions, the underlying 
data can produce MFs with vanishingly small statistical 
error bars, and systematic errors are not reduced by de- 
scribing the data with analytic fits. As a result, we are 
entering the regime where we create more noise than sig- 
nal by comparing analytic descriptions of the MF, which 
impose structure on the data a priori, rather than com- 
paring actual MF datapoints with statistical and system- 
atic error bars. 

Having validated this method to generate a low- 
mass luminosity function from matched SDSS/2MASS 
datasets, future studies (Bochanski et al., in prep.) will 
extend this technique to the entirety of the SDSS foot- 
print, producing simultaneous measurements of the low- 
mass MF and Galactic structure parameters of the thin 
and thick disks. 

7. CONCLUSIONS 

1. We have identified a sample of more than 13,000 
candidate low-mass stars detected within 30 square 
degrees of 2MASS, SDSS and/or GSC imaging. 
Empirical tests verify the sample is more than 99% 
complete to J = 16.2, and more than 90% complete 
at even the reddest colors. 

2. Analysis of the spectroscopic sample reveals that 
exotic contaminants, such as QSOs, extragalactic 
stellar populations, CVs, and white dwarf/M dwarf 
pairs, comprise less than 0.3% of our matched sam- 
ple. 

3. Photometric i — J color provides a reliable proxy for 
stellar spectral type and Teff, we measure disper- 
sions of 1 subclass for a given 0.1 magnitude slice 
in i — J, or of 0.1 magnitude in i — J for objects 
with a given MK spectral type. 

4. In a smaU percentage of cases (20/9624; 0.02%), 
a star's observed spectral type differs significantly 
from that expected on the basis of its photometric 
i — J color. This usually occurs when the 2MASS 
and SDSS detections of a close visual binary are 
improperly matched, and we have identified ani — z 
vs. i — J color cut which can be used to eliminate 
such spurious detections from the sample. 

5. Analysis of the kinematic and spectroscopic prop- 
erties of our sample reveal 88 candidate subdwarfs, 
representing a subdwarf fraction of 0.68% for our 
broader sample. As the bulk of the candidate sub- 
dwarfs (69/88) were identified kinematically, and 
proper motion catalogs are biased towards objects 
with large space motions such as subdwarfs, this is 
an upper limit to the true subdwarf fraction of our 
photometric sample. 
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6. TRILEGAL simulations of the Galactic stellar pop- 
ulation predict giant stars will contaminate our 
sample at the 4% level globally, but with giants ex- 
ceeding 10% of stars with i — J ^ 1.15. Candidate 
giant stars identified in our sample using spectral 
indices verify the colors of giants predicted by the 
TRILEGAL simulation. Restricting our analysis 
to sources with J > 12 and i — J > 1.25, however, 
reduce contamination of our sample by giant stars 
to a managable level (less than 2% in the bluest, 
most contaminated color bin) while sacrificing only 
a handful of bona fide low-mass dwarf stars. 

7. After correcting for unresolved binaries and 
Malmquist effects, we find the luminosity function 
of the Galactic disk is consistent with that mea- 
sured from volume complete samples of the solar 
neighborhood. The dominant systematic uncer- 
tainty in our analysis is due to the remaining uncer- 
tainty in the SDSS/2MASS CMR, which translates 
into a 70% uncertainty at the peak of the stellar 
LF. Systematic effects due to the binary correction 
prescription and the adopted Galactic model are 
comparatively small, at the 10-20% level. 

8. The logarithmically binned mass function mea- 
sured from our sample peaks at 0.3 Mq and agrees 
within la for all but one data point o f the mass 
function measured by iReid et alT ()2002f ) from their 
volume complete sample. A log normal char- 
acterization with A=0.1 (in units of stars (Log 
Mp,)~^pc~ ^ for c omparison with the values re- 
ported by iChabri er 2003), Mc=0.29, and cr=0.32 
provides an adequate description of the data, with 
90% confidence intervals for each parameter of 
A=0.08-0.12 , Mc=0.20-0.50, and cr=0.22-0.37. 
These 90% confidence intervals correspond to lin- 
early binned mass functions peaking between 0.27 
Mq and 0.12 Mq, where the best fit MF turns over 
at 0.17 Mq. a power law fit to the entire mass 
range sampled here, however, returns a best fit of 
a=l.l (where the Salpeter slope is a = 2.35); a bro- 
ken power law returns a=2.04 at masses greater 
than log M = -0.5 (M=0.32 Mq), and a=0.2 at 
lower masses. 

9. We emphasize that comparisions of analytic MF 
fits, rather than the underlying data, can give the 
false impression of MF variations even when the 
underlying data agree well. We are entering the 
regime where we create more noise than signal by 
comparing analytic descriptions of the MF, which 
impose structure on the data a priori, rather than 
comparing the actual MF datapoints themselves. 
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APPENDIX: PHOTOMETRY OF PARALLAX 
STANDARDS ON THE SDSS SYSTEM 



An accurate characterization of the color-magnitude 
diagram for low-mass stars is of critical importance for 
inferring stellar densities and distributions from purely 
photometric datasets. By mining the SDSS database, 
and conducting supplementary observations with the 
NMSU Im telescope, we have obtained colors and mag- 
nitudes on the SDSS photometric system for 76 stars and 
brown dwarfs with measured trigonometric parallaxes. 

Al.l. Mining the SDSS database for Parallax Standards 
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As the SDSS contains photometry of nearly a quar- 
ter of the night sky, thousands of stars with measured 
trigonometric standards he within the survey footprint. 
While the vast majority of parallax standards are so 
bright that they saturate SDSS imaging, parallax stan- 
dards with servicable SDSS observations represent a 
valuable opportunity to map the main sequence color- 
magnitude relation directly from SDSS photometry, obvi- 
ating the need to employ delicate, color-dependent trans- 
formations to place observations from other telescopes 
onto the SDSS photometric system. 

To identify a catalog of parallax standards with 
native SDSS photometry, we have searched the 
SDSS database for detections of objects contained 
within seve ral pu b licly a vailable parallax cat- 
alogs (JMonet et all 119921: Ivan Altena et all 119951: 



Tinnev et al. 1995 : Perr vman et al . 1997: iDahn et aLl 



2002t iGouldfc Chanamel 120041 : fVrba et all 12001^ 
Culling the list to remove objects with flags indicative 
of severe photometric problems (ie, SATURATED- 
CENTER and/or INTERP-CENTER) left 40 objects 
with native SDSS photometry and direct trigonometric 
parallax measurements. SDSS and 2MASS photometry 
for these objects, as well as their measured trigonometric 
parallaxes, are presented in Table El Figure [9] displays 
the locations of these sources in an i — J vs. Mj 
color-magnitude diagram. 

A1.2. Observations of Parallax Standards with the 
NMSU Im 

Additional constraints on the location of the lower 
main sequence in SDSS/2MASS color-magnitude space 
can be derived by transforming observations from sec- 
ondary telescopes onto the native SDSS photometric sys- 
tem. Using time on the NMSU Im telescope made avail- 
able to the public through the NSF PREST program. 



sured trigonometric parallaxes, along with 'deep' obser- 
vations of fields containing late type stars with unsatu- 
rated SDSS photometry. 

The NMSU Im telescope features a 2048x2048 E2V 
CCD providing a plate scale of 0.467 arcsec/pixel. Ob- 
servations were carried out in robotic mode under good 
conditions (clear to light cirrus, seeing 1-2") on Febru- 
ary 6th and 7th, 2007. Images were bias subtracted, flat 
flelded, and overscan corrected before instrumental mag- 
nitudes were measured via aperture photometry, using 
apertures with a 15 pixel radius and sky annuli 15 pixels 
in width with an inner radius of 25 pixels. Obser v ations 
of u'g'r'i'z' photometric standards iSmith et al.l ()2002f ) 
demonstrate internal photometric accuracy of cr=0.04 
magnitudes. 

Comparing instrumental magnitudes to SDSS 
database magnitudes for stars in the 'deep' SDSS fields 
allow the derivation of airmass corrections simultane- 
ously with the zeropoints and color-terms required to 
transform from NMSU instrumental magnitudes onto 
the SDSS photmetric system. The accuracy of the 
NMSU to SDSS transformation is limited by statisti- 
cal uncertainties in NMSU detections of stars in the 
'deep' SDSS fields; residuals in the transformed SDSS 
i magnitudes range from 0.03 to 0.13 magnitudes for 
sources with i = 15 — 18, respectively. Similar results 
are seen for stellar colors, with maximum r — i and 
i — z residuals at the faint (z = 18) end of 0.15 and 
0.2 magnitudes, respectively. As a result, we adopt 
conservative uncertainties for our transformed SDSS 
photometry of 0.15 magnitudes in i and r — i, and 0.2 
magnitudes in i — z. 

The SDSS magnitudes we derive for these sources from 
our transformed NMSU photometry are presented in Ta- 
ble [51 along with their 2MASS magnitudes and measured 
trigonometric parallaxes. These sources are also shown 
in Fig. M 
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TABLE 8 
Native SDSS Photometry of Parallax Standards 



Star 
Name 



SDSS 
i 



SDSS SDSS 



2MASS 
J 



2MASS 
J -H 



2MASS 
H-Ks 



(mas/yr) 



Note. — Reference 1: 'van Altena et all | |1995D 



Reference 2 
Reference 3 
Reference 4 
Reference 5 



Gould & Chanamc i2004D 
Monet ct al L1992.) 
Vrba et al.1 l l2004f 



Dahn et all 1 12002 



ref. 



LHS 1970 


16.26 


0.83 


0.49 


14.58 


0.58 


0.12 


12.9 


0.7 




G 195-022 


12.79 


1.74 


1.00 


10.26 


0.56 


0.30 


64.4 


4.0 




LHS 2364 


17.01 


0.29 


0.16 


15.98 


0.45 


-0.03 


38.3 


2.7 




LP 323-239 


13.55 


1.66 


0.90 


11.18 


0.44 


0.26 


48.1 


1.9 




LHS 2828 


14.81 


2.23 


1.14 


12.00 


0.52 


0.33 


36.5 


0.9 




LHS 2930 


14.29 


2.82 


1.52 


10.79 


0.65 


0.35 


103.8 


1.3 




LHS 3061 


17.54 


1.20 


0.60 


15.52 


0.51 


0.11 


8.9 


0.7 




J090551. 11+553218.5 


13.66 


1.39 


0.78 


11.49 


0.60 


0.23 


20.9 


0.97 


2 


J091841. 52+582747.5 


14.33 


1.44 


0.74 


12.16 


0.58 


0.28 


16.0 


1.2 


2 


J103327.6+472202.6 


14.01 


1.51 


0.83 


11.71 


0.65 


0.26 


15.4 


4.0 


2 


J113212. 96+003632.5 


14.18 


1.44 


0.95 


11.78 


0.52 


0.29 


21.2 


4.0 


2 


J113304.07+131816.9 


14.91 


1.36 


0.69 


12.79 


0.61 


0.27 


10.8 


1.5 


2 


J120723. 99+130213.9 


15.68 


2.01 


1.05 


13.00 


0.60 


0.32 


24.1 


1.4 


2 


J123426. 4+391309. 4 


14.21 


1.36 


0.76 


12.10 


0.57 


0.24 


14.4 


3.1 


2 


J132820.87+300319 


16.73 


2.50 


1.37 


13.31 


0.67 


0.35 


20.7 


2.0 


2 


J135921. 6+251403.1 


15.02 


1.89 


1.00 


12.43 


0.56 


0.31 


24.2 


2.1 


2 


J141436. 71+160121. 7 


14.39 


1.73 


0.90 


12.00 


0.59 


0.31 


28.1 


1.1 


2 


J155719. 67+074500.1 


14.90 


1.44 


0.75 


12.76 


0.52 


0.27 


12.7 


1.5 


2 


LP374-4 


17.01 


0.29 


0.16 


15.98 


0.45 


-0.03 


38.3 


2.8 


3 


LP130-226 


15.32 


2.54 


1.28 


12.20 


0.51 


0.33 


39.6 


1.1 


3 


LP323-239 


13.55 


1.66 


0.90 


11.18 


0.44 


0.26 


49.0 


2.2 


3 


LP499-5 


14.81 


2.23 


1.14 


12.00 


0.52 


0.33 


36.5 


1.0 


3 


LP98-79 


14.29 


2.82 


1.52 


10.79 


0.65 


0.35 


10.38 


1.4 


3 


J003259. 36+141036.4 


23.09 


1.16 


3.69 


16.83 


1.18 


0.70 


30.1 


5.2 


4 


J010752. 46+004156.3 


21.19 


2.02 


2.55 


15.82 


1.31 


0.80 


64.1 


4.5 


4 


J015141. 68+124429.4 


22.85 


1.60 


3.39 


16.57 


0.96 


0.42 


46.7 


3.4 


4 


TVLM 832-10443 


17.08 


2.86 


1.67 


13.13 


0.69 


0.48 


36.0 


0.4 


5 


J074642. 48+200031. 6 


16.08 


2.57 


1.85 


11.76 


0.75 


0.54 


81.9 


0.3 


5 


J082519. 44+211550.3 


20.61 


2.80 


2.76 


15.10 


1.31 


0.76 


93.8 


1.0 


5 


J083008. 16+482847.2 


21.22 


2.35 


3.14 


15.44 


1.10 


0.67 


76.4 


3.4 


4 


TVLM 213-2005 


17.14 


2.77 


1.66 


13.39 


0.65 


0.48 


30.1 


0.4 


5 


LHS 2471 


14.70 


2.60 


1.50 


11.26 


0.60 


0.40 


70.2 


2.1 


5 


J125453. 99-012247.4 


22.25 


2.05 


4.22 


14.89 


0.80 


0.25 


75.7 


2.9 


4 


LHS 2924 


16.19 


2.39 


1.80 
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0.76 


0.48 


95.0 
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5 


LHS 2930 
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2.82 
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0.65 
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16.08 


2.64 


1.79 


11.87 
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0.32 


16.2 


2.6 


4 



Hogg, D. W., Finkbeiner, D. P., Schlegel, D. J., & Gunn, J. E. 

2001, AJ, 122, 2129 
Huchra, J., & Sargent, W. L. W. 1973, ApJ, 186, 433 
Ivezic, Z. et al. 2002, in ASP Conf. Ser. 284: AGN Surveys, 137-+ 
Ivezic, Z. et al. 2004, Astronomische Nachrichten, 325, 583 
Juric, M. et al. 2008, ApJ, 673, 864 

Klessen, R. S., Burkert, A., & Bate, M. R. 1998, ApJ, 501, L205+ 
— . 2002, Science, 295, 82 
Lada, C. J. 2006, ApJ, 640L, 63 
Langon, A., & Wood, P. R. 2000, A&AS, 146, 217 
— . 1992, MNRAS, 256, 641 
Le Borgne, J.-F., et al. 2003, A&A, 402, 433 
Luhman, K. L. 2004, ApJ, 617, 1216 
Lupton, R., Gunn, J. E., Ivezic, Z., Knapp 

Yasuda, N. 2001, in ASP Conf. Ser. 238 

Analysis Software and Systems X, 269 — h 
Lupton, R. H., Gunn, J. E., & Szalay, A. S. 1999, AJ, 118, 1406 



1999, Astronomische 



G. R., Kent, S., & 
Astronomical Data 



Malyuto, V., Oestreicher, M. O., & Schmidt-Kaler, T. 1997, 

MNRAS, 286, 500 
Malyuto, V., & Schmidt-Kaler, 

Nachrichten, 320, 71 
Marshall, D. J., Robin, A. C., Reyle, C., Schultheis, M., Picaud, S. 

2006, A&A, 453, 635 
Miller, G. E., & Scalo, J. M. 1979, ApJS, 41, 513 
Monet, D. G., Dahn, C. C., Vrba, F. J., Harris, H. C., Pier, J. R., 

Luginbuhl, C. B., & Abies, H. D. 1992, AJ, 103, 638 
Morell, O. 1988, in lAU Symp. 132: The Impact of Very High 

S/N Spectroscopy on Stellar Physics, ed. G. Cayrel de Strobel 

&; M. Spite, 519-+ 
Morrison, H. L. et al. 2003, AJ, 125, 2502 
Munn, J. A. et al. 2004, AJ, 127, 3034 
O'Donnell, J. E. 1994, ApJ, 422, 158 
Padoan, P., & Nordlund, A. 2002, ApJ, 576, 870 
Perryman, M. A. C. ot al. 1997, A&A, 323, L49 



SDSS/2MASS M star Lum./Mass Function 



21 



TABLE 9 

NMSU-lM Photometry of Parallax Standards 



Star 


SDSS 


SDSS 


SDSS 


2MASS 


2MASS 


2MASS 


7r 


TT 


TV 


Name i 


r — i 


i — z 


J 


J - H 


H -Ks 


(mas/yr) 


err. 


ref. 


AC +61 23399 9.72 


0.66 


0.32 


8.150 


0.69 


0.14 


32.1 


1.6 


1 


BD+37 2337 9.29 


0.29 


0.12 


8.033 


0.56 


0.11 


24.99 


1.6 


1 


CNLeo 


10.41 


2.39 


1.38 


7.085 


0.60 


0.40 


425.0 


7.0 


2 


G 88-28 


9.66 


0.96 


0.48 


7.818 


0.64 


0.23 


50.9 


2.7 


1 


G96-1 


9.67 


0.16 


0.01 


8.645 


0.46 


0.02 


24.07 


3.1 


1 


G 176-8 


9.76 


0.76 


0.41 


8.107 


0.67 


0.22 


30.1 


2.3 


1 


G 198-37 


9.23 


0.69 


0.33 


7.610 


0.64 


0.20 


43.25 


2.0 


1 


G 199-13 


9.09 


0.22 


0.06 


8.000 


0.50 


0.10 


22.7 


1.2 


1 


G 221-23 


8.37 


0.47 


0.20 


6.873 


0.61 


0.15 


57.53 


2.7 


1 


G 221-24 


9.23 


0.95 


0.45 


7.379 


0.58 


0.23 


57.53 


2.7 


1 


G 222-1 


9.55 


0.71 


0.30 


7.955 


0.72 


0.14 


28.8 


1.4 


1 


G 250-34 


9.67 


0.94 


0.46 


7.872 


0.61 


0.24 


56.51 


2.2 


1 


G 256-10 


10.28 


0.33 


0.07 


9.131 


0.48 


0.13 


27.00 


10.0 


2 


GJ 3452 


9.68 


1.07 


0.54 


7.772 


0.55 


0.27 


76.3 


2.6 


1 


GJ 3504 


15.26 


2.48 


1.40 


12.035 


0.56 


0.34 


50.8 


0.5 


2 


GJ 3855 


14.33 


2.74 


1.50 


10.79 


0.65 


0.35 


103.3 


1.3 


2 


GJ 585.1 


9.15 


0.53 


0.23 


7.693 


0.61 


0.18 


39.0 


1.7 


1 


GJ 9455B 


10.01 


0.19 


0.04 


8.948 


0.38 


0.13 


20.7 


10.6 


1 


GSC 03722-00302 9.70 


0.30 


0.09 


8.450 


0.51 


0.12 


20.39 


2.0 


1 


HD 24916 


7.31 


0.34 


0.12 


6.063 


0.58 


0.15 


63.41 


2.0 


1 


HD 284552 


9.78 


0.38 


0.14 


8.427 


0.59 


0.15 


24.98 


2.0 


1 


HIC 33805 


9.72 


0.64 


0.32 


8.129 


0.59 


0.19 


43.9 


2.2 


1 


LHS 35 


10.68 


1.59 


0.80 


8.424 


0.50 


0.27 


143.0 


4.0 


2 


LHS 197 


14.37 


2.26 


1.14 


11.559 


0.50 


0.29 


51.8 


1.0 


2 


LHS 224 


10.99 


1.75 


0.91 


8.537 


0.44 


0.32 


109.1 


2.8 


2 


LHS 231 
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1.78 


0.90 


11.606 
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0.24 


40.7 
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2 
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0.29 


96.7 


2.3 


2 
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0.87 


0.46 


7.994 


0.60 


0.23 


47.1 


2.4 


1 


LHS 331 


12.17 


1.48 


0.74 


9.983 


0.48 


0.28 


39.9 


1.0 


2 


LHS 2914 
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0.29 
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10.803 


0.47 


0.16 


27.0 


13.0 


2 
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10.35 
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0.61 


0.13 


28.7 


3.4 
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0.32 
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0.70 


0.13 


39.0 
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2.70 


1.71 
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0.64 


0.41 


31.9 


2.9 


3 


Note. — Reference 1: |Porr\ 


anan et al.l (119971) 

( 19951) 














Reference 2: 


van Altena et al. 




Reference 3: 


Tinnev et al.l 119 





Pier, J. R., Munn, J. A., Hindsley, R. B., Hennessy, G. S., Kent, 
S. M., Lupton, R. H., & Ivezic, Z. 2003, AJ, 125, 1559 

Prugniel, P., & Soubiran, C. 2001, A&A, 369, 1048 

Reid, L N., & Gizis, J. E. 1997, AJ, 113, 2246 

Reid, I. N., Brewer, C, Brucato, R. J., McKinley, W. R., Maury, 
A., Mendenhall, D., Mould, J. R., Mueller, J., Neugebauer, G., 
Phinney, J., Sargent, W. L. W., Schombert, J., Thicksten, R. 
1991, PASP, 103, 661 

Reid, L N., Gizis, J. E., & Hawley, S. L. 2002, AJ, 124, 2721 

Reid, L N., Hawley, S. L., & Gizis, J. E. 1995, AJ, 110, 1838 

Reid, N., & Hawley, S. L. 2000, New light on dark stars : red dwarfs, 
low mass stars, brown dwarfs (Springer) 

Reipurth, B., & Clarke, C. 2001, AJ, 122, 432 

Salpeter, E. E. 1955, ApJ, 121, 161 

Scalo, J. M. 1986, Fundamentals of Cosmic Physics, 11, 1 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Schmidt, M. 1968, ApJ, 151, 393 

Sesar, B., et al. 2006, AJ, 131, 2801 

Siegel, M. H., Majewski, S. R., Reid, I. N., & Thompson, L B. 
2002, ApJ, 578, 151 

Silk, J. 1995, ApJ, 438, L41 

Silvestri, N. M. et al. 2006, AJ, 131, 1674 

Skrutskie, M. F. et al. 2006, AJ, 131, 1163 

Smith, J. A. et al. 2002, AJ, 123, 2121 



Space Telescope Science Institute, & Osservatorio Astronomico di 

Torino. 2001, VizieR Online Data Catalog, 1271, 
Stobie, R. S., Ishida, K., & Peacock, J. A. 1989, MNRAS, 238, 709 
Stoughton, C. et al. 2002, AJ, 123, 485 
Tinney, C. G., Reid, L N., Gizis, J., & Mould, J. R. 1995, AJ, 110, 

3014 
Torres-Dodgen, A. V., & Weaver, W. B. 1993, PASP, 105, 693 
Tucker, D. L. et al. 2006, Astronomische Nachrichten, 327, 821 
van Altena, W. F., Lee, J. T., & HofHeit, E. D. 1995, The general 

catalogue of trigonometric [stellar] paralaxes (New Haven, 

CT: Yale University Observatory, — cl995, 4th ed., completely 

revised and enlarged) 
Vanden Berk, D. E. et al. 2005, AJ, 129, 2047 
Vrba, F. J. et al. 2004, AJ, 127, 2948 
West, A. A., Hawley, S. L., Bochanski, J. J., Covey, K. R., Reid, 

I. N., Dhital, S., Hilton, E. J., & Masuda, M. 2008, AJ, 135, 785 
West, A. A. et al. 2004, AJ, 128, 426 
West, A. A., Walkowicz, L. M., & Hawley, S. L. 2005, PASP, 117, 

706 
Williams, C. C. et al. 2002, Bulletin of the American Astronomical 

Society, 34, 1292 
Woolf, V. M., & Wallerstein, G. 2005, MNRAS, 356, 963 
Worthey, G., & Lee, H. . 2006, astro-ph/0604590 
York, D. G. et al. 2000, AJ, 120, 1579 



